ROUGH VOLTERRA EQUATIONS 2: CONVOLUTIONAL 
GENERALIZED INTEGRALS 



AURELIEN DEYA AND SAMY TINDEL 



Abstract. We define and solve Volterra equations driven by an irregular signal, by 
means of a variant of the rough path theory allowing to handle generalized integrals 
weighted by an exponential coefficient. The results are applied to the fractional Brownian 
motion with Hurst coefficient H > 1/3. 
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1. Introduction 



Let x be a general n-dimensional Holder continuous path with Holder exponent 7 > 0, 
an initial condition a G M'^, and cr : M+ X R+ X M"^ ^ R'^''^ a smooth enough function. 
Then a general form of stochastic Volterra equation driven by x (considered as a noisy 
input) can be written as: 



where T an arbitrary positive constant. This kind of system being widely used in the 
physical and biological literature, its noisy version has also been intensively studied when 
the driving motion a; is a Brownian motion [2], [3l [I5] or a general semi-martingale |20) . 
If the coefficient a is also considered as a random function, which is natural in many 
situations, some anticipative stochastic calculus techniques are required in order to solve 
equation and we refer to [H [5], [6l [TH [171 IIS] for the main results in this direction. 
It should be mentioned at this point that the last of those references |T9] is motivated 
by financial models of capital growth rate, which goes beyond the classical physical or 
biological applications of Volterra equations. 

It seems then quite natural to generalize the aforementioned results, and consider sys- 
tems like jl]) driven by general continuous processes, whose prototype can be thought of 
as a n-dimensional fractional Brownian motion. In this case, and when one desires to go 
beyond the Young case 7 > 1/2, rough paths type techniques must come into the picture. 
However, the classical rough path theory introduced by Terry Lyons [I3] (see also the nice 
introductions P, US]) is mostly designed to handle the case of diffusion type equations, 
and there have been an intensive activity during the last couple of years in order to extend 
these semi-pathwise techniques to other systems, such as delay equations |T1| or PDEs 
[H Hi]. The current article fits then into this global project, and we shall see how to 
perturb the original rough path setting in order to handle systems like ([T|). 

Before we come to a description of our main results, let us mention a few choices we 
have made for this paper: 

(i) Like in [T4l \TT\ . we have chosen to work with a variant of the rough path theory 
introduced by Gubinelli in |iO], and called algebraic integration. This method is based 
on some simple enough algebraic considerations, and this relative simplicity makes it 
amenable to intuitions on possible generalizations of the original setting, beyond the 
diffusion case. In the case of Volterra equations handled here, we will see that, in spite of 
the huge amount of technical details involved in our proofs, the main ideas on which our 
constructions rely are quite natural. 

(a) We have specialized equation ([T]) in the following manner: instead of considering a 
general coefficient of the form cr(t, «,?/„), we have assumed that the coefficient a can be 
decomposed under the form (pit — u) cr(x), for a given kernel : ]R_|_ M and a matrix- 
valued function a defined on R'^. Furthermore, an additional hypothesis is made on the 
kernel 0: we assume that it can be written as the Laplace transform of a certain function 
0, namely that 



for a certain /3 > 0. This additional assumption is made in order to take advantage of the 
multiplicative property of the exponential function, and it should be noticed here that the 
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and 
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same kind of results could have been obtained by means of Fourier (instead of Laplace) 
transforms. The integrability hypothesis on 0(^) is morally equivalent to a regularity 
condition on our kernel 0. Once these assumptions are made, and up to a an application 
of Fubini's theorem which can be justified easily in case of a smooth driving process x, 
our Volterra system can be written as: 

yt = a+ di m / e-«(*-") (T(y„) dx^. (3) 

JO JO 

(in) An additional cosmetic change is the following: in order to ease some of our future 
expansions, we transpose the matrix notations given before and set yt = y^. With this 
little change in the usual notations, one is left with the following system: 

/•oo nt 

yt = a+ d^m / e-«(*-")dx„a(yJ. (4) 
Jo Jo 

This is the general form under which we shall solve our Volterra problem. 

With these preliminaries in hand, the main results contained in this paper can be 
roughly summarized as follows (see Theorem 14.161 below for a precise statement) : 

Theorem 1.1. Let x be a n-dimensional fractional Brownian motion with Hurst param- 
eter H > 1/3. Assume that can be decomposed as with [3 = 2, and that a is 
a C^''' -function. Then equation ^ admits a unique solution on any arbitrary interval 
[0,T], in a class of paths called convolutional controlled processes, and where the integral 
with respect to x has to be interpreted as in Proposition ^. 7[ 

Let us now say a few words about the methodology we have adopted in order to solve our 
equation: as mentioned before, it consists in an elaboration of the tools introduced in [TO] . 
Let us recall that these latter reference relies on the definition of an elementary operator 5, 
which transforms for instance a function / of one variable t G [0, T] into a function of two 
variables as {5f)ts = ft — fs- Under some algebraic and analytic conditions, this operator 
6 can be inverted, its inverse is called A, and this inverse allows to construct a generalized 
integral of Young type. If one wants to solve an equation of the form dyt = dxtcr^yt), a 
possible strategy is then the following: remark first that the a priori increments of y can 
be decomposed as: 

{Sy)ts = {Sx)ts Cs + rts, with (s = c^iys), and rts = / dxu [cr{yu) - criVs)] ■ (5) 

J s 

Furthermore, if a; is a 7-H61der process, one also expects y to be 7-H61der continuous. 
Thus, if a is regular enough, ( will inherit the same regularity, and it is also easily 
conceived that r should have the double regularity, namely 27-H61der. This is precisely 
the structure asked in |T0] for the solution to the diffusion-type equation dyt = dxt<7{yt), 
and a process admitting the decomposition (l5l) is called a weakly controlled processes. 
The second important point in the theory is to notice that a reasonable definition of the 
integral dxu [cr{yu) — <^(z/s)] can be given for a controlled process, provided that the so- 
called Levy area associated to x (formally defined as / / dxdx) can be constructed, and 
thanks to the operator A mentioned above. This integration step transforms a weakly 
controlled process into a weakly controlled process, and allows to settle a fixed point 
argument for the resolution of the diffusion equation driven by x. 

Let us try to explain now what has to be changed to the original algebraic integration 
setting in order to handle the case of a Volterra equations: 
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(i) Observe first that in order to solve equation ([4]), the main step is to define accurately 
the rough integral yt{0 = Jq ^^^^^^^^ dxucr{yu) for all t,^ G M+. As will be explained 
at Section 13.11 an important step in this direction is to note that, in order to get some 
increments of y involving only integrals of the form , one has to introduce some twisted 

increments of the form 5yts{0 = ^VtsiO ~ (e~^*^*~*^ — l)ys{0- Then it is easily checked, 
in case of smooth paths y and x, that the a priori twisted increments of the solution y to 
equation (l4|) can be expressed as 

~Syts{0 = j' e'^'^'"'^ dx,a{y,). 

The operator 6 will thus play a central role in our computations, a fact which is reminiscent 
from the calculations contained in [H] for the definition of rough PDEs. It turns out that 
the operator 6 can also be inverted under some algebraic and analytic conditions. This 
inverse gives then raise to a generalized convolutional Young integral, which is at the core 
of our definition of the integral e~^*^*~"^ dx^ (y{yu)- 

(a) The notion of controlled paths has also to be changed for the resolution of equation ([4]), 
and we shall introduce a notion of convolutional controlled path, which will be detailed 
at Section \A7\\ They are basically defined as in equation ([5]), except that 5x is replaced 
by an increment of the form x^, with x\g = (f){t — v)dxy, which is assumed to exist 
once and for all. Then as in the diffusion case, we are able to define a natural extension 
of the notion of integral for those convolutional controlled processes, provided that some 
double iterated integrals based on x can be defined. More specifically, the equivalent of 
the notion of Levy area in our Volterra context is an increment indexed by the Laplace 
variable, of the form: 

xliO= dr]^{r]) e-^^'"''Ux,, e-'^^'^'-^'Ux,,. (6) 

Jo J s J s 

Here again, this definition is only formal in case of a Holder path x, but once it is assumed 
to exist and to satisfy suitable analytic and algebraic hypotheses, a good notion of integral 
can be constructed for convolutional controlled processes. This allows again a fixed point 
procedure in order to solve our Volterra equation. A quick glimpse at the proof of Theorem 
14. 161 will show however that this fixed point procedure is trickier than in the diffusion case. 

(in) An essential step in our approach is thus a good definition of the double integral (l6|), 
and the study of its regularity in This can be done quite easily (up to some Garsia 

type regularity theorems which have to be proven) when x is taken as a Brownian motion, 
and when the integrals with respect to x are interpreted in the ltd sense. However, an 
important part of the current article will be devoted to the definition of ([6]) when a; is a 
fractional Brownian motion with !/?> < H < 1/2. Here again, the path we have followed 
in order to obtain this definition is not completely standard, and let us say a few words 
about it. Indeed, the usual way to define a double integral like ^ in case of a fBm 
is to use Stratonovich integrals, in the sense of the Malliavin calculus as explained in 
|16| . However, this way to compute our iterated integrals involves a decomposition of the 
Stratonovich integral into a Skorokhod type term plus a trace term, which is hard (though 
not impossible) to analyze in case of an exponentially weighted integral like ours. We have 
thus decided to adopt another strategy, and have resorted to an analytical approximation 
of the fractional Brownian motion introduced in [23]. This latter approximation, which 
will be recalled at Section [3 has the advantage to yield almost explicit and elementary 
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computations, based on the analysis of singularities for some locally analytic functions 
defined on the complex plane. In our opinion, the calculations we obtain are thus more 
elegant than in the Malliavin calculus setting. 

As we mentioned before, the resolution of rough Volterra equations relies thus on a 
few simple ideas. These ideas are however long to formalize when one wishes to give 
most of the details of the calculations, which explains the bulk of the current article. It 
should also be mentioned at this point that we could have tried to solve equation ([I]), 
in its general form, without recurring to twisted convolutional increments as we did, by 
just following the standard algebraic integration formalism. This is in fact what we have 
done in the companion paper [7J, and this idea works fine for the Young case, namely for 
a Holder coefficient 7 > 1/2. However, this first method fails to give a global existence 
result in the case 1/3 < 7 < 1/2, mainly because the Picard iterations don't lead to 
a contraction property (we refer the reader to [7] for a more detailed argument). This 
important drawback justifies the introduction of the convolutional generalized integration 
we have used in the current paper. 

Here is how our article is organized: we recall some basic definitions of algebraic inte- 
gration at Section [2l Section [3] is devoted to the simpler case of Young equations, which 
allows to explain our method with less technical apparatus. Then at Section [4] we move 
to the rough case of our Volterra equation, and explain all the details of the method we 
have chosen in order to solve it. Finally, we apply our theory to the fractional Brownian 
motion case at Section [5l 

2. Algebraic integration 

This section is devoted to recall the very basic elements of the algebraic integration 
theory introduced in ^0], in order to fix notations for the remainder of the paper. We 
also include a proof of the existence of the so-called sewing map A which is simpler than 
the one contained in the original paper [10], and is even a further simplification of the 
proof proposed in |11| . 

2.1. Increments. As mentioned in the introduction, the extended integral we deal with 
is based on the notion of increment, together with an elementary operator 6 acting on 
them. The notion of increment can be introduced in the following way: for two arbitrary 
real numbers £2 > ^1 > 0, a vector space V, and an integer /c > 1, we denote by Ck{V) 
the set of continuous functions g : [£1,^2]^ V such that gti- tk = whenever ti = tj+i 
for some i < k — 1. Such a function will be called a (A; — 1) -increment, and we will set 
C^{V) = Uk>iCk(y). The operator 6 alluded to above can be seen as an operator acting 
on /c-increments, and is defined as follows on CkiV): 

k+l 

6 : C,{V) ^ C,^,{V) {5g)t,..,,,, = Yy^9t,-u-t,,,, (7) 

i=l 

where tj means that this particular argument is omitted. Then a fundamental property 
of (5, which is easily verified, is that 55 = 0, where 55 is considered as an operator from 
Ck{V) to Ck+2{V). We will denote ZCk{V) = Ck{V) n Ker5 and BCk{V) = Ck{V) n lm5. 

Some simple examples of actions of 5, which will be the ones we will really use through- 
out the paper, are obtained by letting g & Ci and h G C2. Then, for any t,u,s G [^1,^2], 
we have 

{S9)ts = gt- gs, and {5h)tus = hts - htu - Ks- (8) 



6 



AURELIEN DEYA AND SAMY TINDEL 




Furthermore, it is readily checked that the complex (C*,5) is acyclic, i.e. ZCk+i{V) = 
BCk{V) for any k > 1. In particular, the following basic property, which we label for 
further use, holds true: 

Lemma 2.1. Let k > 1 and h G ZCk+i{V). Then there exists a (non unique) f G CkiV) 
such that h = 6f . 

Observe that Lemma [2TT] implies that all the elements h G C2{V) such that Sh = can be 
written as h = 6f for some (non unique) / G Ci{V). Thus we get a heuristic interpretation 
of S\c2{v)- it measures how much a given 1-increment is far from being an exact increment 
of a function (i.e. a finite difference). 

Notice that our future discussions will mainly rely on fc-increments with k < 2, for 
which we will use some analytical assumptions. Namely, we measure the size of these 
increments by Holder norms defined in the following way: for / G C2(y) let 

sup iM|_, and C^{V) = {feC2{Vy,\\f\\,<oo}. 
In the same way, for h G CsCV), set 

\\ht \\v 

||/i|Lp = sup -. ^ (9) 

P,,fi-P,] h = ^hi,0<p,<fi 

i 

where the last infimum is taken over all sequences {hi G C^^V)} such that h = hi and 
for all choices of the numbers pi G (0,z). Then ||-||^ is easily seen to be a norm on C^iy), 
and we set 

C^iV) := {heCsiVy, \\h\\^<oo}. 

Eventually, let Cl~^(y) = U^>iC3(V^), and remark that the same kind of norms can be 
considered on the spaces ZC^iV), leading to the definition of some spaces ZC^iV) and 
ZCl'^{V). In order to avoid ambiguities, we shall denote by N[f] Cj] the K-Holder norm 
on the space Cj, for j = 1, 2, 3. For ( G Cj{V), we also set M[C] Cj{V)] = sup^g[^^.£2]j ||Cs||y- 

With these notations in mind, the following proposition is a basic result which is at the 
core of our approach to path-wise integration: 

Theorem 2.2 (The sewing map). Let p > 1. For any h G ZCi^{[0, 1]; V), there exists a 
unique Ah G ([0, 1]; \^) such that 6 (Ah) = h. Furthermore, 

\\Ah\\^<c^^^[h;c^{v)], (10) 

with Cfj_ = 2+2'^ Xlfcli This gives rise to a linear continuous map A : ([0, 1]; V^) — »■ 
([0' 1]; V) such that 6A = Idzc!^([o,i]-v)- 

Proof. The original proof of the existence of A (with a somewhat different constant C^) 
can be found in [10] and has been simplified in |11) . We give here a more elementary 
proof, which will be easy to adapt to the pertubated incremental operator 6 (see Section 
Q. 

Uniqueness. Let M, M' G such that 6M = 6M' = h. In particular, 6{M - M') = 0, 
and thus, according to l2.1l M — M' = 6q, with q G Ci. But then g G Cf with > 1, hence 
q is constant, and as a consequence, M = M' . 
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Existence. By definition of ZC^^, we know that there exists B G C2 such that SB = h. 
Consider now the sequence (7r"')„ of dyadic partitions of [0, 1], that is 

tt" - {0 = < < • • • < = 1}, with = ^, 
and set, for all s,t & [0, 1], 

if tt" n {s,t) = 0, 
Ks ={Bts- Btjs - if TT" n {s, t) = {t^}, 

Bts - bI^, - Bun - Y!r}^ B^^^ if TT" n {s, t) = {t]<---< tf}. 

It is readily checked that the mapping M" : s,t Mj" is continuous on [0,1]^. We 
are now going to show that the sequence {M"-)nen converges in the space C{[0, 1]^; V) of 
continuous functions on [0, 1]^, endowed with the norm J\f[- ; ^^(V^)]. 

Let then s, t G [0, 1], n G N, and denote 

7r"n(5,t) = {t;<t;Vi<---<^r} 

= {t^;'<i?i<---<^St2<^r'}, withj</<2«. 
If s < and t < then 

7r"+^ n (s, t) = < < < • • • < i^-i < t^h 

and in that case 

i=j 

which, since 5B — h, leads to 

We proceed likewise for the cases (s > t2j'\,t > t^^^), {s < t2j'\,t > t^^^) and (s > 
iJ_\,^<^S+\),to finally get 



n 

+ 



21^ 1X2/^/ V 2^- 

and thus 

Since we have considered /x > 1, this proves that the series ^„A/'[M"''"^ — M"; C2(^)] 
converges, and thus Xln(^^"^^ ^ ^'^") converges in C([0, 1]^; endowed with the norm 
A/'[- ; (^2(1^)], the latter space being complete. But, invoking the fact that M° = 0, we 
have = En=o (^"^^ - ^"), which entails the uniform convergence of towards 
an element M G C([0, 1]^; V). We can already notice that for all n, MJ^ = 0, which yields 
the same property for M, so that M G C2. 
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Take now 0<s<M<t<land denote vr" n {s,u) = {t], . . . tt" n [u,t) 
{t}, . . . , tp}, hence tt" n (s, t) = {t], ...,tf}U {t^„ tf,}. Thus, 



i-i 



I'-i 



— Bf o — Bfn — Bftn — \ Btn 4-n — N Bfn 4-n — B+n, xn . 
tS •!> tl;/ / J ^i+l^i / J ^i+l^i ^jl^l 

i=j i=j' 

We will assume that t", > u, the case t"; = u leading to the same relation (fTTll . Then 



MZ = Bts + 



Bus — Bt^s — Butv- — 7 ^ Btn t" 

J I / ^ z+1 i 



+ 



I'-i 



Btu — Bt^u — Btt"^ — > Bt^ 



-Bwts + ~ Bus — Btu — B, 



3 ' 



which can be written as 



Ml = M:, + M^^ + htus-ht^,ut-. 



(11) 



Since h G C3, lim„_^oo ^t",utp = 0, so that, by letting n tend to infinity in the previous 
relation, we get 5M = h. 

Finally, let us show that for any s,t E [0, 1] and n G N, 

\\MZ\\v<c^Af[h;Ci;iV)]\t-sr, (12) 

which will prove inequality (fTOj) and as a consequence, the Holder regularity of M. To this 
end, fix s,t G [0, 1], n G N. If vr" n (s,t) = 0, the result is obvious. If vr" n (s,t) = {t]}, 
Ml = i5B)ujs = hu^s, hence ||M,"J|v < Af[h; C^iV)] \t - s]" < c,U[h; Ct^iy)] \t - s\r If 
vr" n (s, t) = {t], . . . , tp}, pick A; G {j + 1, 1} such that 

l^fe+i " ^k-i\ ^ I _ j _ I 1^ ~ -^1 • 

At this point, the previous relation does not seem very relevant insofar as the distances 
between two successive points of tt" are equal. In fact, this relation will make sense when 
we iterate the scheme. Consider indeed the new partition it = {t", . . . , t^+i, . . . , tf} 

and define M/^ according to the same principle as M/J, using n instead of vr" fl (s, t). Then 



ikf-f" — M^It — Bfn fn — BfTifn — Bfn fn — flfn fUfn 



and as a result 



\Ml-Ml\\v<mh;Ci;{V)]- 



2M 



We iterate the procedure until the partition reduces to the empty set, to get 

WM^Wv < m-, c^iv)] \t-sr (^2+2^ < c,M[h- c^,{v)] \t - s\ 



□ 



The following corollary gives a first relation between the structures we have just intro- 
duced and generalized integrals, in the sense that it connects the operators 5 and A with 
Riemann sums. 
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Corollary 2.3 (Integration of small increments). For any 1-increment g G C2{V), such 
that 6g e set 6f = {Id-M)g. Then 



j=0 

where the limit is over any partition Tits = {to = t-, ■ ■ ■ i^n = s} of [t, s] whose mesh tends 
to zero. The 1-increment 5 f is the indefinite integral of the 1-increment g . 

Proof. For any partition Ilf = {s = to < < ••• <'tn = t} of write 

n n n 

{5f)ts = '^{Sf)u+iu = ^9u+ik -^K+iui^9)- 

i=Q i=0 i=0 

Observe now that for some fi > 1 such that 6g E Cl^, 

n n 

II - < ^imy, c^iv)] \t-s\, 

4=0 i=Q 

and as a consequence, liming, |^o X]r=o (^5') = 0. □ 

2.2. Computations in C*. For sake of simplicity let us assume for the moment that 
V" = M, and set Cfc(M) = Ck. Then the complex (C*, 6) is an (associative, non-commutative) 
graded algebra once endowed with the following product: for g E Cn and h E Cm let 
gh G Cn+m the element defined by 

(5'^)tl,...,im+n-l = 5'il,...,in^tn,.--,im+n-l' ^1 5 • • • 5 ^ [^l,^2]■ (13) 

In this context, we have the following useful properties. 

Proposition 2.4. The following differentiation rules hold true: 

(1) Let g^h be two elements ofCi. Then 

5{gh) = 6gh + g6h. (14) 

(2) Let g E Ci and h E C2. Then 

6{gh) = 6g h + g 6h, 5{hg) = 6h g — h6g. 
Proof. We will just prove (fT4ll . the other relations being equally trivial: if g, h E Ci, then 

[^i9h)]t, = gtht - gshs = gt ih - hs) + {gt - gs) h, = gt {Sh)^^ + {6g)^^ h,, 
which proves our claim. 

□ 



The iterated integrals of smooth functions on [^1,^2] are obviously particular cases of 
elements of C which will be of interest for us, and let us recall some basic rules for these 
objects: consider f,gE C^, where is the set of smooth functions from [£1,^2] to R. 
Then the integral / dg f, which will be denoted by J{dg f), can be considered as an 
element of . That is, for s,t G [^1,^2], we set 

Jts{dg f) (^j dgf^ = dgju- 
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The multiple integrals can also be defined in the following way: given a smooth element 

h e and s,t E [ii,i2\j we set 



Jts{dgh) = 




In particular, the double integral Jts{,dpdp f^) is defined, for G Cf^, as 

Jtsidfdf f')={^j dfdf = 1^ dfl {df f) . 

Now, suppose that the nth order iterated integral of df"' ■ ■ ■ dp still denoted by J{df" 
■■■dp p), has been defined for p,p ...,p e C^. Then, if G C^, we set 

Jtsidn'dp ■ ■ ■ dpp) f dp+' X, {dp ---dp p) , (15) 

J s 

which defines the iterated integrals of smooth functions recursively. Observe that a nth 
order integral J{dp ■ ■ ■dpdp) could be defined along the same lines. 

The following relations between multiple integrals and the operator 5 will also be useful 
in the remainder of the paper (see e.g. [TTj for a proof of these elementary facts): 

Proposition 2.5. Let f,g be two elements of . Then, recalling the convention ( fJ^j . 
it holds that 

6f = Jidp, 6 iJidgf)) = 0, 5 iJidgdf)) = i5g)i6P = J{dg)J{dp, 
and, in general, 

n-l 

6 {j{dP ■ ■ ■ dp)) =Y^j{dp^^^ dp+') J {dp ■ ■ ■ dp) . 

i=l 

3. VOLTERRA EQUATIONS IN THE YOUNG SETTING 

Recall that we wish to solve equation ([41), and we start this program by studying the 
Young case, i.e. the case of a driving process x which is assumed to be 7-H61der continuous 
with 7 > 1/2. This allows us to introduce most of the general tools used in the sequel, and 
this section is thus conceived as an introduction to the rough case which will be treated 
later on. In order to get a feeling of the kind of structure needed in order to deal with 
Volterra equations, we will start with some heuristic considerations, which are basically 
justified in case of a smooth driving noise x. Then we shall proceed to define rigorously 
the equation, and solve it in a suitable class of functions. 

3.1. Heuristic considerations. Assume for the moment that x is a smooth process, in 
which case equation ([4]) is well defined and solvable when a is a regular coefficient. In 
order to get an intuition of the natural operators associated to our equation, let us recast 
it, in quite a redundant way, as a system: 

yt = a + VtiO d^ , . 

m =/o*e-«(*-^)rfx.a(i/.). ^ ^ 

Notice that the first relation above is not sufficient in order to determine as a function 
of y. However, the second one defines y without ambiguity. 

As we already mentioned in the introduction, if one wants to generalize the system 
we have just written to a non smooth signal x, it is now easily seen that the main step 
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is to give a rigorous meaning to the integral e~'^^^^'"^dx,u a{yy) defining ytiO- To this 
purpose, and having introduced the main tools of algebraic integration in the last section, 
the first idea one may have in mind is to get a suitable expression of the increments 
{6y{^))ts = y{Ot ~ y{Os for a given ^. And indeed in case of a smooth driving process x, 
invoking equation ( flGl l. those increments can be written as: 

mO)ts = f e'^^'~^Ux, a{y,) + ats{0 f e-^^'-'^^dx, a{y,) 

J s Jo 

\~^^'--'"^dx,a{y,) + ats{i)Ui), 

where we have set ats{C) = e"^^*""*^ — 1. Notice now that the first term e~^^^~'"^dxv(j{y.u) 
above is really similar to what one obtains in the diffusion case, namely an integral of the 
form J^. However, the second term dtsiOVsiO ^ little clumsy for further expansions. 
Hence, a straightforward idea is to make it disappear by just setting {Sy)ts{C,) = {Sy)ts{0~ 
ats{^)ysiO- Then the last equation can be read as {5y)tsiO = Jg ^^^^^^^^''dx^ a{y^), and 
the system ([T6l) becomes 



yt =a + ytiO 0(0 d^ , . 

with the initial condition yo = 0. This very simple fact, together with the nice algebraic 
properties which will be seen below, converts the elementary operator S into the central 
object in order to solve our Volterra system. 

These preliminaries being admitted, we shall essentially focus in the sequel on the 
process y, by merging the two equations of the last system into a single one: 

mts{0= I e-«(*-^)dx,afa+ / dr^^{r^)Uri)] ■ (18) 



The original solution process y can then be recovered in an obvious way, and we shall 
solve the Volterra equation under the form (fTSl) . 



3.2. Convolutional increments. Let us turn now to the main concern of this section, 
that is the definition of a complex which behaves nicely for the definition of our 

Volterra problem. 

Notice that, due to the fact that e"'^'-*^"*^-* is nicely bounded only for ti > t2, our 
integration domains will be of the form Sn = iS„([£i, £2]), where iS„ stands for the n- 
simplex 

Sn = {(tl, . . . , tn) : 4 > tl > t2 > ■ ■ ■ > t„ > h}. 

Let then be a separable Banach space. In order to define the basic family of continuous 
increments we will work with, we first need to specify the (Banach) functional space each 
yts{-) will belong to, with a special emphasis on the Laplace coordinate. In fact, the 
calculations to come (see for example Lemma [3^ incite us to consider the £^-type space 
induced by the norm 



Jo 
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where /3 > is fixed. Then, we define Cn,^ as the space of continuous applications from 
Sn to CjsiV). Observe that an operator 6 : Cn,f^iV) Cn+i^piV) can be defined just like 
in (I7l). In particular, if A G Ci^/^iV) and B G C2^/3(y), the relation (l8|) is still valid. 

As we have seen at Section 13.11 a suitable family of operators related to our Volterra 
equation is given by 6 : Cn,p{V) Cn+i,f3{V), defined for any positive ^ by 

= - for A G C^AV), (19) 

where G Sn+i- In the remainder of the paper, we will explicitly write the 

variable ^ down only when there might be a confusion. Thus, we simply write 6A = 
SA — aA, where we made use of the convention (fT3l) . As in Section [2Tn one can define, 
for n > 1, 

ZCnAV) = CnAV) n ker(5), and BCnAV) = CnAV) n Im(5) 

In fact, when = R'^ or = M'^''^, endowed with their natural Euclidian norms, the 
convention (fTSll can be extended according to the following principle: 

Lemma 3.1. Let M G C„,/3(K^'') and L G C„(R'). Then ML, defined by the relation 

belongs to Cm+n-i,i3{^^) ■ Moreover, when n = 2, the following algebraic relations hold 
true: 

S{ML) = SML-MSL, and 6{ML) = 6ML-M6L. 
Proof. The first part of our claim is an obvious consequence of 

As for the algebraic relations, the first one follows from Proposition 12.41 while 

= [{^M)t^t^t^ - at,t.,Mt2t,]Lt^...t^^, - {M6L)t^,„t^^^. 

□ 

With these preliminaries in hand, it is now easily shown that the perturbed operators 
6 preserve some important properties of the original coboundary 6: 

Proposition 3.2. 66 = 0. More precisely, the couple (C*^^(V^),5) satisfies ZCnA^) ~ 
BCnAy) for aUn>0. 

Proof. This proof is borrowed from [H], Proposition 3.1] and is included here for sake of 
completeness. If F G Cn,f^(y), according to the fact that 66 = and thanks to Lemma 
13.11 we have 

66F = {6-a)[{6-a)F] = 66F-6{aF)-a6F + aaF 
= -6a F + a6F-a6F + aaF = aaF-6aF. 
Furthermore, it is readily checked that 

{6a)tus = atu CLus, {t, u, s) G S3, 

which gives 66 F = 0. 



ROUGH VOLTERRA EQUATIONS 



13 



The fact that Im5|j^^^ ^^y) = ker ^(y-, can be proved along the same lines as for 

the (C*,5) complex [10]: pick A G Cn+i,i3(y) such that 6A = 0, and set Bt^,,,tn = Ai...t„s, 
with s = 0. Then 

Thus, setting C = we get 5C = A. 

□ 

The cochain complex {C^:^i3{V), 6) will be the structure at the base of all the constructions 
in this paper. Let us also mention at this point that, when the meaning is obvious, we will 
transpose the notations of Section [2TT] to our convolutional setting. Furthermore, whenever 
this doesn't lead to an ambiguous situation, we will write Cn,p instead oiCn,/3iV). 

We will now define an equivalent of the iterated integrals of Section 12.21 in our convo- 
lution context: for two smooth functions f,g,ii<s<t<i2 and ^ > 0, define 

Jts{dgfm = j^e-^^'-'^Ug^U 



and for h e CP° 



2 ) 

t 



Jts{dghm= I e-^^'-^'Ug^Ks. 

J s 

Once these elementary blocks have been defined, the iterated integrals 

J{dg\..dg-f) for g\ . . . , g^J e , (20) 

should be defined as functions of several variables, according to the same recursive prin- 
ciple as in Section 12. 2t 

JUdg' . . . dg- . . . , D = e-^^'^'-'^Ugl J^dg^ . . . dg^ • • • , D- 

In particular, J'ts{dxdx){^,r]) = e'^^^~'"^dxy e~'^^"~'^^dx^. 

The following relations between 6 and these integrals will be useful for our purposes: 
Proposition 3.3. Let f,g e . Then 

S [j{dg /)) = ~5 [j{dg 5/)) = J{dg) 5f. 
Proof. Straightforward. 

□ 

3.3. Holder spaces and A-map. In the Young setting, it will be enough to let our 
solution live in some Holder-type spaces. Indeed, one expects the solution y to (fTTl) to 
belong to a space of the form Cf (as defined at Section [2TT1) for any 1 — 7 < /? < 7, where 
7 is the Holder regularity exponent of the noise x. Since /? + 7 > 1 in this case, the 
exponentially weighted integrals with respect to x can be interpreted in the Young sense, 
as will be explained below. 

As far as the path y alluded to in (fTTl) is concerned, we also expect his increments 6y 
to be regular enough. Thus, we shall resort to the following natural Holder spaces: 

^2,/3 ■= {y ^ ^2,/3 : ^f[y; C2J ■■= sup ■^J^*^'^^^] < 



0<s<t<T 



t-s 
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:= {y E C,,p : ~5y G C>^}. 

Notice that our definition of the space Cf^ is based on the twisted operator 5 instead of 
5. For any h G C^^p, set, just as in the standard case, 



M[h-C^^f]:= sup 

0<s<u<t<T \t - U\' \U - S]"^ 



M[h- := inf 1^ Ar[/i,; C^^f-''\ /i = ^ /i„ < < /i| , 

where the last infimum is taken over all sequences {h E C^,^} such that h = J2i for 
all choices of the numbers pi G (0, /x). Denote also ZC^ ^ := lm{6) HC^^ and observe that 
the property ZC2 = {0} if p > 1 remains true for 6: 

Lemma 3.4. If p>l, then ZCi^^^ = {0}. 

Proof. Let M = 6f e ZC^j^. Consider the telescopic sum {Sf)ts = T.to ^'~''^'K^ fk+iU 
with respect to the partition Uts = {s = to < ti < ... < t^ = t} of the interval [s,t]. Then 

n 



i=0 



which tends to as the mesh \Ilts\ of the partition decreases to 0. 

□ 

As we already mentioned, an essential tool in order to define generalized convolutional 
integrals is the following inverse of the operator 6: 

Proposition 3.5 (The convolutional sewing map). Let /i > > 0. For any h G ZOi^ p, 
there exists a unique Ah G Cl^^ such that 6{Ah) = h. Furthermore, 

Af[AkC^J<c,Af[kC!;j, (21) 



with = 2 + 2^^ YlT=i W- This gives rise to a linear continuous map A : ZCi^n — >■ C,^ 



2,/3 



such that 5 A = Id^ri^ ■ 

^'-3,13 



Proof. It follows the same line as the proof of Proposition I2.2[ 

Uniqueness. Let M, M' G Cl^j^ such that SM = 6M' = h. In particular, 6{M - M') = 0, 
hence M - M' e ZCl^j^. Thanks to lemma [33, we deduce M = M'. 

Existence. As in the standard case, consider B E C2 3 such that 6B = h and construct 



if 7r"n (s,t) = 

{SB)u^s ■dn^nis,t) = {t^} 



and if vr" n (s, t) = {t] < ... < 

MZ = Bts - Buy - E ^-'^'^K.U - e ^'"^^ B^ 



S5 

J 
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where stands for the n-dyadic partition of [0, 1]. 

It is readily checked that M" is continuous from [0, 1]^ to Cf^. Moreover, if for instance 
7r« n = {t^+i < < - < ^2/"^^}, with s < t^+\ and t > then 

i-i 

i=j 

and since 6B = h, this yields: 

«'-M«;r,.^{.(l)\(J,)"}. 

This estimation remains true for the other cases of intersection of vr" with Using 
the same arguments as with A, we thus get the existence of a limit M of M" in €2^/3. 

The fact that 6M = h can be proved just as in Proposition 12.21 and it is the same for 
the estimation 

□ 

We also have the following equivalent of Corollary [231 which links A with convolutional 
Riemann sums, at our disposal: 

Corollary 3.6. For any 1-increment g G €2^/3 such that 6g G Cg^ (l^ > set 6f = 
{Id-M)g. Then 

n 

{6f\ = , lim e-'^~''^'^9u^.u in Cp, 

1=0 

where the limit is over any partition Uts = {to = t, . . . ,tn = s} of [t, s] whose mesh tends 
to zero. 

Proof. We use the same arguments as in the standard case, starting from {Sf)ts = 
E^e-^'-'^^^KSfK.u- 

□ 

3.4. Young convolution integral. Recall that, according to the notations of Section 
13.21 the Volterra equation ( fTSl) we are interested in can be read as 

yo = 0, 6y = J (^dx ^ + _^ dV kvMv)^ ) • (22) 

We will now define integrals of the form J{dxz)^ such as the one appearing in the right 
hand side of equation (l22l) . when x, z are only 7-H61der with 7 > 1/2. This will rely on 
the following assumption, which is trivially met when x is a smooth path: 

Hypothesis 1. Assume that, for some 7 G (1/2, 1), x is a path in C1{E}''^), allowing to 
define an increment x^ G C2 ,y(ffi^''^) which satisfies 5x^ = 0. 

Remark 3.7. The increment x^ represents morally the integral J'{dx), which will be defined 
as a Wiener integral in the fractional Brownian case (see Section [3^ . Furthermore, under 
Hypothesis [H the increment xj^ = dC, is well defined as an element of C2 . 

The fact that x^(-) G will be simply ensured by the condition /q°°(1 + ^''^)|0(OI dC, < 00. 
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Theorem 3.8. Let x be a path from [0, T] to M}''^ satisfying Hypothesis [H for a given 
7 e (1/2, 1). Let z G CI, and for ^ G R+, set 

J{dx z){^) = S:\^) z + A[xUz]{^) = {Id-A6)[x^ z]{^). (23) 

Then 

(1) J'{dx z) is well defined as an element ofC^.^, and coincides with the usual Riemann 
integral /* e~^^^~'"^ dx-u Zy when 7 = 1. 

(2) For a constant Cx > 0, we have, for all £i < £2, 

^^[Jidxz);Cl^i[ium < Cx {Ar[^;C°([£i,£2])] +eW[^;C7([£i,£2])]} , 

where the norms Af have been defined at Section \3.3l J\f[z;C^] := sup^^<g<;^2 Iksll 
and e = — £i\- 

(3) It holds that, for any ii < s < t < £2, 

n 

Jts(dxz)= lim y^e"-^*"*'+i^x! . Zt. in C^, 

1=0 

where the limit is over all partitions Hts = {Iq = t, . . . ,tn = s} of [s, t] as the mesh 
of the partition goes to zero. 

Proof (1) In the regular case, the integral J'ts{dx z)(^) = e^^^^"""^ dx^ Zy is defined in 
the Riemann sense, and it is readily checked that 

Jtsidx z) = x]^ Zs + Jts{dx 6z), (24) 

and hence 

Jts{dx 6z) = Jts{dx z) - x]^ Zs. 

Applying 5 to both sides of this last relation and taking into account Proposition 13.31 
Lemma [3TT] and Hypothesis [U we obtain: 

5 {^J{dx5z)^ = —Sx^ z + x^ 5z = x^ 5z. 

Now, if and z are 7-H61der continuous with 7 > 1/2, A can be applied to the relation 
above, and one can write Jts{dx5z) = A{x^6z). Plugging this equality into (l24ll . we 
obtain the expression (l23ll . Thus our integral coincides with the usual one in case of a 
regular process x. 

Since 27 > 1, the item (2) is a direct consequence of the contraction property (l2Tl) of 
A. As for (3), it stems from Corollary 13.61 

□ 

3.5. Volterra equations. We are now ready to solve equation (1221 ). by interpreting the 
integral J'{dx a{a + drj (j){ri)y{ri))) as in Theorem I3.8[ Before stating the main theorem 
in this direction, let us introduce the subspace C°'^ of Cj.^ induced by the norm 

where A/'[^; C°^] := suPq^^^j^ Af[ys', C^]- With this new space in hand, we can prove the 
following elementary lemma, which will be used throughout the proof of the theorem: 
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Lemma 3.9. Let o E Cl and for any y E C°'^, set y = a + dr](j){r])y{7]). Then 
a{y) G Cj and 

J\fHy)-Cl]<c^Af[mC',::i]. 
Moreover, if y^^\y^'^^ G C°'^([£i, £2]) are such that yf^ = yf^, then 

M[a{y^'^)-a{y^'^)-Cl{[l,A,])]<c^[l+M[^^^^ 
with e = \ii — £2\- 

Proof. The three inequalities are mostly due to the obvious estimation 



Indeed, we have for instance 

POO 

^f[6icriy)Uv < ikiloo/ dmowmtsioi 

Jo 

^ POO POO 

< Ik'lloo \t-s\^ [ATlm C7,,] + Ar[y; C° 

and therefore Af[a{y);Cj] < \\a'\\^Af[y;Cl'^]. 

The second inequality can be obtained in the same way, after noticing that, for all 

SG [£1,^2], 



7 n 



>7J 



Jo 

As far as the third inequality is concerned, we can invoke the classical estimation 

- cr{y^'%s\\v < IWlUSiy^'^ - y^%s\\v + \W"\\oomy^%s\\v 

{\\yl''~y?'\\v + \\yi''~yi'%)- 

□ 

We are now in position to prove the 

Theorem 3.10. Assume HypothesisUl holds true for some 7 > 1/2, and that a G C^'''. 
Then equation admits a unique solution in C°'^, where the integral J'{dxa{a + 

drj (l){r])y{ri))) stands for the Young convolutional integral introduced in Theorem \3.8[ 

Proof. Let e > (we shall fix this constant retrospectively), / G N, and suppose that we 
have already constructed a solution y^ G C°'^([0, Ze]). If Z = 0, then y^ = jjo = 0. We 
mean to extend into a solution y'''^^ G C°'^([0, (/ + l)e]), by resorting to a fixed point 
argument. 

Step 1: Existence of invariant halls. Let y G C°'^([0, (l + l)e\) such that y\[o,ie] = y^ and set 
z = T{y) the element of Ci,^([0, + defined by 5|[o,Z£] = y'' and for all s,t E [0, {l + l)e], 
{Sz)ts = Jts (^dxo-{y)^^ where, as in Lemma IXOl y = a + dr](j){r])y{7]). 
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We know from Theorem 13.81 that 

M[i-A,,{[le. {I + 1)^])] < c. {Ar[a(y);C?([0, (/ + l)e])]+e<U[a{y)-Cl{[^, {I + l)e])]] , 
which, together with Lemma [3^ gives 

m Cl^ille, {I + 1)5])] < {l + C°':^([0, (/ + l)e])]} . 

li < s < le < t < {I + l)e, we use the relation 55 = to deduce 

and hence 

< 2max(Ar[5;C7,,([/e,(/ + l)e]),Ar[y';C7,,([0,fe])]) \t - . 

Furthermore, for all s,t E [0, (/ + l)e], Zg = {Sz)so, and thus 

Af[~^;Cl{[0, (/ + 1).])] < Ar[5;C°'^^([0, (/ + 1)5])]T^. 
We are therefore incited to set 

Ni+i = max (2(1 + nAr[y';C7,([0,fe])], 44,^(1 

Indeed, for such values, it is readily checked that if N'[y]Ci'l^{[0, [l + l)^])] < Ni+i, 
then X[z;Cl^{[0, {I + 1)e])] < ^ and ^^[z;Cl^{[0, (/ + 1)^])] < which gives 

^^[z;Cll^{[0, {I + 1)£])] < Ni+i. In other words, the ball 

= {y e cT,v^([o, {I + i)e]) : mie] = y\ m-A?M {i + 1)^])] < iv^+i} 

is left invariant by V. 

The independance of e with respect to the initial condition allows to repeat the 
scheme with the same e and thus to get a sequence of radii {Nk)k>i such that the sets 
Q^k).^ are invariant by T. Of course, the definition of the latter mapping has to be adapted 
(in the natural way) to each of those sets. 

Step 2: Contraction property. We will now search for a division of the previous intervals 
[/£,(/ + l)e] into subintervals \^e,le + rj\,[le + rj^le + 2rj\, . . . of the same lenght (possibly 
depending on e, /), on which a contraction relation holds. 

For I G {1,2}, let y« G CT;:^([0, fe + r/]) such that yfl,^^ = y\ A^fyW; C°'^([0, fe + 
^])] < and denote zS""^ = r(y*^*)), where F is defined as in Step 1, but restricted to 

C°'^([0, le + f]]). According to Theorem 13.81 

Ar[5« -z(2);C7,,([k,fe + r/]) < c,,4Ar[a(y«) -a(y(2));C?([fe,k + r/])] 

+ r^W[aiyW)-aiy^%Cji[le,le + v])]}, 

which, together with Lemma [3.91 implies 

Af[~z^'^ - C7^([k, le + v])] < {1 + iV,+i} r;W[y(i) - y^^); Cl'^^^Hle, le + r,])]. 
Since the processes y^^^ — y^'^\ z^^^ — z^"^^ vanish on [0, le], we can more simply write 
Ar[5(i) - C7^([0, le + v])] < {1 + iV,+i} v'^^'^ - y^'^; C7,,([0, le + v])]- 
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Besides, {z^^^ - z^^)), = - z^^)),,;,, so that N[i^^^ - i^^); co^([o, fe + r^])] < M[i'^^^ - 

5(2); Cj^([0, le + r7])]r/^. Finally, we get 

Ar[5(i) - 5(2); C°',^([0, fe + r,])]< cl^ {1 + iV,+,} (1 + T^)r^W[y(^) - y^^); C;%^([0, k + r^])]. 

Fix then rj = inf (e, (2c^ {1 + A^;+i} (1 + T'^))^^/'^) . In this case, F becomes a strict 
contraction on the set 

{y e C°;,^([0, le + v]) : moM = y\ CTv'([0, le + r^])] < iV,+i}. 

Using the stability of Qyl^^i^^^y, we can easily show that the latter set is invariant by F too 
(cf Lemma [3.111 below) . Consequently, there exists a unique fixed point in this set, which 
we denote by y''^. Since 77 does not depend on y', the same calculation then remains true 
on the (invariant) set 

{y e C°',^([0, le + 2r/]) : y|[o,ze+,] = t\ + 2r^])] < iV,+i}. 

Thus, y^'^ can be extended in a solution y''^'' defined on [0, le + 2rf\ and proceeding so until 
the whole interval [/e, (Z + l)e\ is covered, we get the expected extension 

□ 

Lemma 3.11. With the notations of the preceding proof, the set 

{y e cl;{[o, le + v]): mis] = y\ cTv'([0' + ^])] < ^'+1} 

is invariant by F. 

Proof. Let y an element of this set and z = T{y). Set 

^ _ Ivt ift <le + r] 

^* ~ \e-(*-('^+^))yi,+, if te[le + r], {I + l)e]. 

Then y is easily shown to be continuous, that is y G Ci,y([0, (/ + l)e]). Moreover, if 
s, te[le + ri, (/ + l)e], {5y)ts = 0, whereas if s < k + < t, {5y)ts = e-^''^^'+'^^\5y)ie+r,,s, 
so that Af[y;Cl^{[0, {I + l)e])] < Ar[y;Cl^{[0,le + v])]. Since Ar[mCl,{[0, (/ + 1)^])] < 
X[y; C;%([0, le + r;])], we deduce Af[y; ^"'.^([O, (/ + 1)^])] < Af[y; C?;,^([0, le + v])] < iV,+i, 
which means that y G Q^'l/^^)^- But we know from the first step of the preceding proof 
that Q^'J+i)^ is invariant by F, and so, if z = F(y), Af[z;C°{^'^{[0, (/ + l)^])] < A^'^+i. It is 

now clear that z = z\[o,i£+n], which finally leads to A/'[z; C°'^([0, le + r]])] < J^lz; C°'^([0, [1 + 
l)e])]<Ni+,. 

□ 

3.6. Application to fBm. We now aim at proving that the previous results can be 
applied to a fractional brownian motion X = (X^^^, . . . , X^")) with Hurst parameter 
H > 1/2. Before we start with this program, let us recall what we mean by fBm in this 
paper (we refer to [H] for further details on this process): for computational sake for the 
case 1/3 < H < 1/2, we will consider X as a centered Gaussian process indexed by M 
(even if our equation is indexed by [0,T]), with covariance 

Rnit, s),, = ^(X«xi^')) = i5,,(|s|2^ + |tp^ -\t- s|2^), s,teR. 

We assume that the underlying probability space {Q, JF, P) on which X is defined is such 
that fl is the Banach space of all the continuous funtions Co(M; R"), which vanish at time 
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0, endowed with the supremum norm on compact sets. P is the only probability measure 
such that the canonical process {Xf, t G M} is a n-dimensional fBm with parameter H 
and the cr-algebra is the completion of the Borel cr-algebra of Q with respect to P. 

In order to apply our general results to the fBm, we need to define Wiener integrals 
with respect to X. To this purpose, denote by H the completion of the M'^-valued step 
functions £ with respect to the inner product 

n 

((l[o,ii], • • • , l[o,t„]), (l[o,si], • • • , l[o,s„])) = ^RHisi,ti), Si,ti e M. 

i=l 

When H > 1/2, it can be checked that this inner product can be expressed as: 

{fi9)n = CHf2[ [ fi'^gi'Uu-vr^'dudv, with CH = H{2H-1). (25) 
Jr Jr 

for all f,gEH. It can then be shown that the family of Wiener integrals {X{h)] h G Ti.} 
with respect to X forms an isonormal Gaussian process, with E[X{hi) X{h2)] = {hi, h2)n- 

With these notations and facts in hand, a natural definition of X^ is as follows: for 

^, s, t G M+ and i < n, set 



XlfiO= fe~^'-'-''^dXj^^:=X{h{t,s-0), with K{t, s; = e-^'^'-''h[,,y 

J s 



V Ci 



(26) 

where Cj denotes the i^^ vector of the canonical basis in M". We will now show that 
this process satisfies Hypothesis [H under some integrability assumptions on 0. First, it 
is readily checked that 6X^ = 0. In order to prove that X^ G Cj,^ a.s, we shall use a 
Garsia-Rodemich-Rumsey (GRR in short in the sequel) type result, which is an extension 
of the original paper |9] in 3 directions: (i) Like in [10], we will get a regularity result for 
a general function R defined on ^2, which is not necessarily the increment of a function 
/ G Ci. (ii) The conditions on R involve SR instead of SR (iii) R also depends on the 
Laplace variable ^. It should be noticed at this point that, for the remainder of the section, 
1S2 stands for S2{[0,T]). 

Proposition 3.12. Let (V, ||.||) a Banach space and fix ^ > 0. Let _R : ^2 x V such 

that R.XO ^ C2(V") and define 



0<v<w<T 



(l){\w-v\) 



where ip,(t> : are strictly increasing functions and 0(0) = 0. Assume now that 

there exists some C(^) > such that, for all s < t in [0,T], 

sup \\mum<r' f 1 ^it-s). (2?) 

s<u<t \\t — S\ j 

Then, for allO<s<t< T, 
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Proof. See Appendix. 



□ 



A preliminary step, before we prove the desired continuity result for X^, is to show 
that this process is at least an element of C2{M}''^) for any fixed ^. This is achieved in the 
following lemma: 

Lemma 3.13. The process defined by formula admits a modification X^'* such 
that, almost surely, X]'*{^) G C2([0,T]) for any ^ > 0. 

Proof. Let us lean on the following version of the Kolmogorov continuity criterion: con- 
sider a process {Zts{^), s,t e [0,T], ^ e M"*"} living in any of the Wiener chaos associated 
to X, and assume that for all X G N and all (52,^2,6) ^ [0,^]^ x [0,N], 



E 



- Zt,s,m\' < CN {\si - s,r + - + lei - ^rn , (28) 



for some 01,02,03 > 0. Then Z admits a modication Z* such that, almost surely, Z*[^) 
is continuous for any ^ > 0. 

To show that X^ actually satisfies (i28ll . suppose for instance S2 < si < t2 < ^i- Then 



SI 



ti 



t2 



S2 



But, on the one hand, relation (12511 yields 



E 



g-«i(ti-«)rfXJ|2 



t2 



and likewise 



E 



ch 



Sl 



S2 



t2 Jt2 



< C Iti - t2\ 



,2H 



< C Isi - S2 



\2H 



On the other hand, still evoking relation (i25l l. we have 



E 



/"'[e-6(ii-«) _e-«2(*2-«)] dXu\ 

J Sl 



<c I ' du f dv |e-«i(*i-^) - e-«2(*2-")| 



\u — V 



' Sl ^ Sl 

which, together with the estimation 

^Q-Clih-u) _ g-6(t2-«)| < |g-Cl {*!-«) _ g-Cl(t2-M)| _|_ |g-6(*2-«) _ g-6(i2-«)| 

< ^l\tl-t2\ + \t2-u\\^^-^2\ 

< iV|ti-t2| + T|ei-6l, 



leads to 



E 



Jsi 

We are thus in position to apply our general results to the fBm case: 



□ 
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Theorem 3.14. Let X he a fBm with Hurst parameter H > 1/2 and 7 G (1/2, if) such 
that /o°° rf^ 10(01(1 + D < 00. Then the process defined by the Wiener integral 
satisfies HypothesisU\for 'y a.s. Consequently, if in addition, a G C^'*, the system 



' {5Y\s = Jts [dX a (a + d^ 0(O>^(O 



1^0 = 

admits a unique solution in Cj''^ a.s. 

Proof. As mentioned before, we just need to check that G C'^^^^iM}''^) . Further- 
more, with Lemma 13.131 in hand, and thanks to the fact that 6X^ = 0, we can apply 
Proposition [3H to X\ with ^(x) = x^p and 0(x) = x'^+^/p, to obtain ||X4(0II < 
c\t — s\'^ (f^7,2p(0)^^^^5 where p is an arbitrary strictly positive number and 

Therefore, we just have to prove that A/'[f/^\'5'; C^] = t^^0(O(l + D(^7,2p(O)^^^^ < 00 
a.s, since in this case, X^ G C2,-y and N'[X^]C2^^] < cA/'[f7^''2jf ! '^7] < ^ ^-S- fact, we 
are going to show that E[Af[u!^^^2p'i 



< 00. 



Let us start with the Jensen inequality 

POO 

Jo 

Notice then that, as we already mentioned in the proof of Lemma 13.131 

E[\\xuorn<c\w-vr\ (30) 

which leads to E[(j^^2piO] — ^ frp^^rplw — vf^^~'^'^'^~^'^ dwdv. This means that if we take 

7 G {1/2, H) and p > 1/{H — 7), E[U^^2p{0] — ^ foi^ some constant M independent of 
^, and as a consequence, 

noo 

E[M[U][^^; < / di 10(01(1 + e). 

Jo 

The proof is now easily finished. 



□ 



4. Rough Volterra equations 



Our aim is still to solve equation ^Tf} or (l22l) . in a case where x satisfies Hypothesis [H 
but where we replace the condition 7 > 1/2 by 7 > 1/3 only. Like in the Young case, 
our first task is thus to give a suitable interpretation to the integral in (l22ll . which goes 
beyond the Young case. We will see that the key to this improvement is to introduce a 
new class of processes. 

4.1. Convolutional controlled paths. As in the Young case, let us start with some 
heuristic considerations: let us go back for a moment to the Volterra equation under the 
form 

yt = a+ / (f){t - v) dxya{y^), 
Jo 
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and assume that x is a smooth path in Cf^ and a, are regular coefficients. Then the 
equation above admits a unique solution y, whose increments can be decomposed into: 

{^y)ts = (t){t- v) dxy a{y^) + [ - v) - (j){t - v)] dxy a{yy) = x]^ a{ys) + ns, 
with 







t 

V 

t 



rts = I (j){t-v)dx^{5a{y))^s+ I [(t)it - v) - (j){t - v)] dx^ a{y^) := r]^ + r'^ 



ts- 



This is exactly the structure which will be imposed for the solution to our equation, and 
let us analyze it a little further: if we assume now that x has only a regularity of the form 
Cj with 7 > 1/3 and that y is for any k < 7, then we expect x^ to be an element of C2 , 
under some regularity conditions on (which will in fact be assumed to be a differentiable 
kernel). As far as the remainder term r is concerned, we expect it to inherit the Holder 
regularity of y and x for r^, and the regularity of for r^. Hence, the remainder term 
r should be an element of C|''. It is also worth recalling from Remark [3.71 that, if x is a 
path allowing to apply Fubini's theorem, then x^ should satisfy: 



f 00 

^1= I xliO md^- (31) 
Jo 

It is thus natural to formulate the following assumption on our driving process x: 



Hypothesis 2. Assume that, for some 7 G (1/3,1/2) and j3 > 0, x allows to define a 



process x^ G C2 «(M^'"') such that 6x^ = 0. Set then 



00 



xuomd^- 







Remark 4.1. Notice that, contrary to Hypothesis [H the index (3 and the exponent 7 may 
be different here. In fact, for some computational reasons that will arise in the proof of 
Theorem 14.161 we shall be prompted to take (3=1. Therefore, from now on, let us only 
focus on the spaces C7i(M^'") (i G {1,2,3}), that we more simply denote by C7(R^'"'). 

Fix an interval I = [a,b] C [0,T] and denote e = \I\ = b — a. With the above 
considerations in mind, the natural spaces to work with in order to solve equation (fTTl l 
can be defined as follows: 

Definition 4.2. A path y G Ci(/;M'^) is said to be a convolutional process controlled by 
x^ (with regularity k) if 6y can be decomposed into: 

{6y)ts = [xl CsY + rts, with ( G C^{I; M"''^), and r G Cl^{I; R'). (32) 

Denote the space of such controlled paths by Q^{r, M'^) and for any h G M.'', write Q1{I; M'^) 
= {y e Q(/;]R^) : ya = h} . Then the norm associated to Q'^{I;R'') is 

Af[y; Q^]=M[y;C^] +Af[C,C'i] +^f[C■,C^] +^f[r■,Ci% 
Notice that 2/ 1/3 < k < 7, Q'^ C Cj and 

Ar[y;C^{I)]<c,e^-^Ar[y; Q^{I)]. 
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In fact, as in the Young case, we shall focus on the form (fTSl) of the original equation 
which also involves a process y indexed by the Laplace variable ^. In this setting, the 
same reasoning as above applied to y leads to the introduction of the following spaces: 

Definition 4.3. A path y G Ci{r,M.^) is said to be a process controlled by (with 
regularity k) if 6y can be decomposed into: 

i5y)ts = (xl CsY + hs, with C e Cr(/; M"''^), and r G Cl'^I'. 
where we recall that we have set C\ := C\y according to Remark [^Tlj Denote this second 
space of controlled paths by Q'^{I]M!') and for any h G Ci, write Q^(/;R^) = {y G 
Qr{I-W) ■ ya = h}. Then the norm associated to Q''(J;]R ) is 



2ki 
2 J 



It is then readily shown that the space of controlled processes is stable by composition 
with a smooth enough function: 

Proposition 4.4. Let z G Q''(/;M'=) with decomposition (Elj, a G C2'^(M'=;M') and set 
z = c^z). Then z G Q'^{T,M}), and it can be decomposed into 

6z = (x^O* +f, 

with 

(s = C{Da{zs)y, hs = Da{zs) rts + [6{o{z))ts - Da{zs){6z)ts\ , 

where Da stands for the matrix-valued coefficient (|^)i<i</,i<j<fc, and the norm of z as 
a convolutional controlled process can be bounded as: 

U[z- Q^il; R')] < {1 + ^^[z; M'^)]^} . (33) 

Furthermore, if z^^\z^'^^ G Q'^{I) are such that z^a^ = z^a\ then 

Af[a{z'^'^) - criz^% Q\I- M})] < c,,,(.),,(.) M[z^'^ - z^'^- Q^{I; R% (34) 

where 

< {1 + Ar[;.«; Q^i^R")] +Af[z^'^; Q''{I;R%' . 
Proof. It is exactly the same as the proof of [TOl Proposition 4], replacing 6X with x^. 

□ 

Finally, let us mention that, in the remainder of the article, we will write Q'^{I) and 
Q^(/) instead of Q'^(/; V) and Q^(/; V) whenever this does not lead to an ambiguous 
situation. 

4.2. Integration of controlled processes. We now aim at giving a precise sense to the 
integral Jts{dxa{y)) which appears in (l22l l and stands for e~'''^~'^^ dxuu{yu) in case of 
smooth processes. As 7 < 1/2, we can no longer resort to Young's interpretation. In fact, 
in order to define this integral, we will rely, as usual in the rough path theory, on the a 
priori existence of some Levy area type process adapted to our problem (notice that the 
following hypothesis covers Hypothesis El): 

Hypothesis 3. Assume that, for some 7 G (1/3, 1/2), x allows to define three processes 
x^ G C7(Mi'"), G C2^(M"'") and x^ G C3'^(M"'") satisfying ~6x^ = and 
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Remark 4.5. In case of a smooth process x, the increment represents now the double 
iterated integral e'^^^'""^ dx^ x^^,, which can also be written, with a slight adaptation 
of the notations of Section 13.21 as the (partially) integrated Levy area 



^tsiO = dr](/){r])Jts{dx (g) dx){^,r]). 
Jo 

As for x^, it is given in this case by 
^tusiO = / e-^^^-'^Ux^ (g) {6x^)yus = / e'^^^^''^ dx,o ® / [0(^; - w) - (f){u - w)]dx^ 



The last ingredient we need before we can integrate convolutional controlled processes 
with respect to the increment dx is a matrix equivalent of Lemma [3?T1 ii A,B G M^'', 
denote A- B = Tr{AB*). Obviously, \A ■ B\ < and hence: 

Lemma 4.6. If M e C2(M'''0 and L e C„(M'='0, then M ■ Le C„+i(M) and 

6{M ■ L) =6M ■ L - M ■ 6L. 

Here is now the natural way to integrate convolutional controlled processes in our 
context: 

Proposition 4.7. For two given coefficients 7, k such that 1/3 < /t < 7, let x be a process 
satisfying Hypothesis{^ Furthermore, let z G Q'^(/;M") with decomposition 

{5z)ts = {xlQ* + where ( G C^{F, M"'"), r G C^,^{F, M"). (35) 

Define A by Aq = h (where h & Ci) and 

{6A)ts = x^z + x^-C + A(5V + ■ {6Cy - ■ O- (36) 

Finally, set 

J{dxz) = 6 A. 

Then: 

(1) A is well-defined as an element of Q''(/;]R), and Jts{dx z){S,) coincides with the 
integral e~^^^^'"^dxv z^ in case of two smooth functions x and z. 

(2) The semi-norm of A in Q'^(/;M) can be estimated as 

H\A- M)] < c.,{{M\z; Clil- M")] + M")]} , (37) 

for a positive constant depending only on x. 

(3) It holds 

n 

Jts{dx z) = \im [s^i+uU + ^l+uU ■ Q] «^ ^1' (38) 

for any ii < s < t < £2, where the limit is over all partitions Hts = {to = 
t, . . . ,tn = s} of [s, t] as the mesh of the partition goes to zero. 

Remark 4.8. It is certainly possible to state and prove continuity results for our extended 
integral in terms of a sequence converging to x in the sense of convolutional controlled 
processes. We did not go into these considerations for sake of conciseness. 
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Proof of Proposition \4. 7[ z,x are two smooth functions, then J'{dx z) can be defined 

as a Riemann integral, and as in the Young case, one can write: 

Jts{dxzm= [ e-^^'-'^Ux, z, = £\0 + f e-^'^'""^dx,{5z),s. 



Plugging the decomposition (i35l l for {5z)ys into this last expression, and observing that, 
thanks to some elementary matrix manipulations, we have dx-u (x^^ Cs)* = {dx-u ®xlg) ■ (*, 
we end up with: 

JUdxzm = S:\Ozs+ [ e-^'^'-^Ux^lixlCr + r.s] 



J s 



or otherwise stated: 



Jts{dx r) = Jts{dx z) - x\^Zs - xl^ ■ Q- (39) 

In order to analyze the term Jts{dx r), let us apply, like in the Young case, 5 to both mem- 
bers of the equality above. This gives, owing to Proposition [331 Lemma [46l Hypothesis^ 
and using the decomposition ([35l) : 

~5{j{dxr)^ = xUz - 6x^ ■ C + S:^ ■ {SO* 

= x^r + x^-{6Cy-x^-C 

When all these terms have a Holder regularity greater than 1, we are now in a shape to 
apply the operator A, which gives: 

J{dxr) =A{x^r + x'^- {60* - x^ ■ C) ■ 

Plugging this equality back into ( l39l ). we have proved the relation 

J{dx z) = £'^ z + x"^ -C + ^{x^r + x'^ ■ {60* - ■ C) , 

in case of some regular functions x and z. 

(2) Let us analyze the two terms of the remainder R oi A defined by (l36l) . namely: 

R = x'^C + Hx^ r + £2 . ^SQ* - £^ ■ C*)- 
For the first term, we have 

M[i' ■ c-A^{i)] < e'^^-^^M[i'-A\i)]m-A{i)]- 

As for the second term, we use the contraction property (l2Tl) of A to deduce: 

j\r[K{x'r + x'-{6o*-x^-cyA''{i)] 

< £W[A(£^ r + £2 . (SO* - x^ ■ Cy, C2"^''(/)] < Ce^ {1 + 11 + HI) , 

with 

I = f^[x\ CT+^(/)] < ^[x'; C]{I)] ATfr; (/)], 

II = Af[x' ■ {60*;Cl^^\I)] < e^-''Af[x';Cl\I)]Af[C,C^{I)], 

III = Af[x^ ■ C; C^^^^I)] < e^''^'-''W[x^; C^^I)] ATfC; C°(/)]. 
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Thus, we get Af[R;C^'^{I)] < e'^^^^'"'^ jV[z; Q'^{I)]. Besides, we already mentioned that 
N'[z;C^{I)] < c^e^-'' Af[z; Q^{I)]. The estimation ([37j) easily follows. 

(3) Remark that 

SA= {ld-AS){xh + x^ ■(*)■ 
The result then stems from Corollary I3.6[ 

□ 

In order to define the term J'{dxa{y)) in our Volterra equation, we need the following 
multidimensional version of the previous proposition: 

Definition 4.9. We say that z E Ci(]R'^'') is controlled by (with regularity n) if z = 
{z'^^\ . . . , z^'-^), with z^''^ G Q^(R''). Denote Q'^(R''''-) this set of controlled processes, and 
define, for any z G Q''(]R^''), 

\i=i 

Let us also introduce the set Q''(M*^'^) along the same principle, together with the norm 

\i=i 

Corollary 4.10. If z e Q'^(M"'''), the process J{dxz) (with values in M!) defined by 

J{dxzY'^ = J{dxz^^) for allie {1,...,/}, 

belongs to Q^iM}). Moreover, the conclusions of Proposition^^ still hold in this context. 
In particular, 

Ar[J{dxzy,Q^{[i,,i2\;R')] 

< {(A/'[z;C°([£i,£2];ffi"'')] + \i2-iir''f^[z; Q'^([£i,^2];M"'')]} • 

4.3. Localized controlled processes. In order to get a global solution for our rough 
Volterra system, we still have to perform a technical step. Indeed, like in the Young 
case, we will solve the equation by patching solutions defined on small intervals, and 
this patching procedure will involve a localization of some convolutional paths around a 
certain smooth increment /, which represents in general an initial condition. The current 
section is thus devoted to adapt our previous definitions and propositions to this localized 
setting. Notice that we assume, throughout the section, that x satisfies Hypothesis [3l 

Fix thus an interval / = [a, b] and denote e = b — a. The following subsets of M.^) 
will come into play: 

Definition 4.11. Let f G C2(/;M''). A process y G C7(/;M^") is said to be K-weakly 
controlled around f if 

- fts = (xLO* + rl, with C e C^{I; M"''^) and r-y G C^'^lJ; M^). (40) 

Denote A'j^^{I]Ml') the set of K-weakly controlled around f processes such that Za = h, 
and for any y G ^j^(/;M'^), define its semi-norm by: 

M[y-,Al,{I-,R')]=^[Cy;C',{I)]+Ar[Cy;C^An]+^[^^ 
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The following elementary facts are worth noticing: obviously, ^q^(/) = the 
norm A^[.; A^i^il; R'')] coincides with7V[.; Q^{r, and for any / G C^il; M^), Al,^{I) C 
Q'^{I). The important point in our localization around / is precisely that this latter 
increment does not play any role in the computation of A^[?/; M'^)] (thus the new 

notation A^, instead of A/", for the norm of y). 

Let us now see how the spaces Aj f^{I) pop out naturally when one integrates a convo- 
lutional controlled process. 

Proposition 4.12. Let y E Q^(/;R'') with decomposition Sy = (x^C^)* + f^. Set y = 
ao + /o°° dUiO m)-^ Then y G Al^{I; M^), with fts = d^m atsi^e^^^'-^HO and 
h = ao + (/)(^) h{^). Moreover, 

M[y;Al,{I)] < c. {Afly; Q''{I)] + e'-^Afik C,]} . (41) 
Proof. U s < t E I, write 

POD 

i5y)ts = I d^mmtsio 

dUiomtsio+ / d^m^^tsiouo 



poo poo 

= (4C|)*+/ da{OrUO+ d^matsmy)saiO + fts. 
Jo Jo 

Set a = CI, rl = /,°^ dUiO {rliO + ^tsiOmsaiO}- Then 

AA[r^;Cr] < c{Ar[P;Cn +Ar[mCn} < cAr[y; Q% 

and \\{6y),,\\ < \\ U + \t ~ M[x';QWCy + \t - sf'' Af[ry;Cn But ||^|| < 
\t- s\ Af[h; £1], hence Af[y; Cf] < e^'''Af[h; Ci] + c^Af[y; Q^], and (HI]) is thus proved. 

□ 

An analog of Proposition 14.41 concerning the composition of a localized controlled pro- 
cess with a smooth function is the following: 

Proposition 4.13. Let y E A'f h{I), and consider a function a E C^'^ . Then a{y) E 
-^Daih)/ a{h)(^) ^'^^ httvc thc following bound on the norm of a{y): 

(.)(/)] 

< c.,. { 1 + Ai [y; Al^{I)f + e'-^M [y; Al^{I)]M[f; Cl{I)] + Cl (/)] } . (42) 

Moreover, ify^'^\y'-^^ E Ajj^{I), 

A-[a(y«)-a(y(2));Q«(/)] 

< c.,.Ar[y(y - Q«(/)]{1 + M[y^'^;Al,iI)r + M[/'^; ^^,,(/)]^ 

+ s'^-Aflf; Clmi + A'[y(i); Cr(/)] + A-[y(2); c^^I)])}. (43) 



ROUGH VOLTERRA EQUATIONS 



29 



Proof. If s, t G /, write 

S{o-{y))ts = / dXDa{ys + X{Sy)ts){Sy)ts 
Jo 

= [ dXDa{y., + X{6y)ts)fts+ [ dX Da{y, + X{5y)ts){{xlQr + rl) 
Jo Jo 

= /^^W^ + (xLC(^))* + r^)•^+r,^^'^ (44) 
with a^'^ = QDa{ys)\ 

= f dX [Da{ys + X{5y)ts) - Da{ys)] {x^QT + f dX Da{y. + X{6y)ts)rl, 



r^)''= / dX[Daiys + Xi6y)ts)-Daiys)]fts + [Daiys)-Daiya)]fts. 
Jo 

By standard computations, 

ATlC^y^; C°] + f^lC'-'^; Cf] + ^[r-^y^''; C^] < c.,. {l + M[y; Al,]'} . 

Besides, 

llr^^-'ll < \\D'a\\oo^f[y;C^]^^[f;C',]{\t-s\'^^ + \s-ant-s\}, 

and hence A/'[r'^(J')'^; Cf*^] < CaJ^[y; A'j^f^jN'lf; C^je^''^ ■ Finally, going back to decomposi- 
tion fl44l). we obtain: 



mcjiy))ts\\ 

< \\Da\U \t - s\ Af[f; Cl] + |t - s\'' U[x\ Cl]M[y; ^^,,] \\Da\U + \t- sj''^ AT^^^); C^], 

sothat A/'[a(?/);CJ'] < c,,,^ {e^-'^TVl/; Ca^] + [y; + 7V[r'^(«'); C|'^] }, which achieves the 
proof of (Il2|). 

As for (l43l) . we have, with the notations (l44ll . 
If we refer now to the proof of [10, Proposition 4], we effortlessly get 

As far as r^^y^^^^'"^ — r'^^*'*^')'^ is concerned, notice that 
\\rt'''^'-rt'''^'\\<\t-s\M\!-A^)^ 

dX \\Da{y^}^ + X{by^^\,) - Da{y^}^) - Da{yf^ + X{by^\,) + )|| 

+ \\Da{y^^) - Daiy^:^) - Z}a(yf ) + Dai/- 



Some standard computations (see e.g. [2T|, Lemma 3.1] for further details), using differ- 
entiations along the path 
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defined for ^,1^ E [0, 1], then lead to 

< c^s'--^^[f■,Cm {1 +Ar[yW;Cr(/)] +Ar[y(2);Cr(/)]}Ar[y« -y(2).c^^(J)]. 
Inequality (j43ll easily follows. 

□ 

Observe again that A'Jj^ is a subset of Q'^, which means that, for any path z G A'j f^, 
the integral J'{dx z) is defined thanks to Proposition I4.7[ In the particular context of a 
process z G A.'j f^, the bounds on this generalized integral can be improved as follows: 

Proposition 4.14. If z E ^j_/j(/;M"), then the semi-norm of the process z in Q''(/;M) 
defined by Za = h E Ci and 6z = J{dxz) can be estimated as 

U[z- Q^{I)] < c,, {U\z- C?(/)] + e^-''M[z- Aly,{I)] + Cl{I)]] . (45) 

Proof. According to Proposition 14.71 z can be decomposed as a controlled process, with 
= z and P = f^'^ + P'^, where 

p'i = x2.C^ and P'2 = A(£V + x2. 5(^-^3. + 

First, since {5z)ts = fts + x]X! + rfs, 

^[C;C^iI)] = Ar[^;Cf (/)] < {e'-^W;Cl{I)] + 6''-^M[z;Al,iI)]} . 

As for the remainder term, we have A/'[r^'^; Cf"] < c^e'^^'^^'^^ Ai[z; Aj^i-Xl)]^ while, thanks 
to the contraction property (fTOj) . 

Afr''; Cr(/)] < c. {e^M[z; Al^il)] + CUl)]} . 

Finally, Af[{SS)ts; C^] < {\t - M[z; ^^,,(/)] + \t - s\'+^Af[f; Cl{I)]}, hence 

Ar[5;Cr(/)] < {e'^-^M[z;Al,{I)]+e'+^--^[f;Cl{I)]} , 

which achieves the proof of (l45ll . □ 

Remark 4.15. If / G C2(/;M'''') and G M'''^ we can define ^^^^(1; M'''') along the same 
lines as Q'^(J; M''''). If 2; G ^^^(1; M"''), inequality jlSl) remains true, that is 

Af[2; Q%I;R')] < {ATf^; C°(/; R"'')] + ^^''^A^l^; ^^,(/; M"'0] + ^'"''W; C2i(/; M"-')]} , 

(46) 

where z is defined analogously to Proposition 14.141 

4.4. Rough Volterra equations. We are now in position to prove the main result of 
this section: 

Theorem 4.16. Let 7 G (1/3, 1/2) and 1/3 < k < 7. Assume x satisfies Hypothesis IM 
and a G C3'''(Ri''^; M"-'^). Then Equation ^ admits a unique solution in Qg([0, T]; R'^). 

Proof. As in the Young case, the solution we are looking for is seen as a fixed point of 
some naturally defined application F. The fixed point argument is then divided into two 
steps: we first establish the invariance of some well-chosen balls of Ql, and then show a 
contraction property on these balls. 



ROUGH VOLTERRA EQUATIONS 



31 



Step 1: invariant balls. Fix a positive integer and consider a sequence of intervals 

so that [0,T] is covered by a 



rN 



[C , C+i] with C = and = £ 



N 



]N _ ]N 

''n 



finite union of (J„ )„>o. Introduce also a sequence of balls 

Bt = {ye QVn ) ■■ = hn, Cf. = ct{K), f^[y; Q'^iC)] <{N + nD 



n] 



We are thus given a control over 



T{y) is the path in C]{I^' 



where hn G £i is such that Af[hn] Ci] < {N + 
both 6y and the initial condition jjiN. If y E -B^", 
defined by the two conditions: Zin = hn and for all s,t E , {Sz)ts = Jts{dx aiy)), with 
y = aQ + d^4>{C)y{C)- With these notations, and using the previous propositions, we 
are going to prove the existence of two constants q;i,Q!2 > such that the sets B^" are 
invariant by F and the following property holds: 

If y e 5^", then M[yiN_^^\Ci] < {N + n + l^K 



H 



Thanks to (H), the local solutions can then be patched together, as we shall see at the 
end of the proof. 

Let y e bI^, ~z = F(y). As in Proposition [Ha denote fl", = d^4){0ats{0e'^^'~''^ 
hn{Cj and /i^ = ao + order to estimate N[z] use successively 

(il), (iD and (IlID, together with the fact that 

to get 



+ e 



1+7-2K 



N[y; Q«(Zf )]Ar[/i„; C,] + el~^N[K; C,]]. (47) 



If one desires to stay in the ball i?^" after applying F, one is naturally led to consider 
the system 

2a2 - (7 - k) < ^2 
2ai - (7 - K + 2(1 - k)) < a2 
ai + ^2 - (1 + 7 - 2k) < 0:2 
ai — (1 — k) < 0^2, 



(48) 



which reduces to 



02 < 7 — K 

«1 — 1 < ^2 ~ 't- 



In fact, for some reasons that will arise soon, we should add the conditions 02 < ^-y^ and 
«! — 1 > ^2 — 7, which turn the previous system into 

< a2 < 2^ 

^2 — 7 < «! — 1 < tt2 — 

Notice that the conditions above can be easily met (and are assumed to be met in the 
sequel) whenever k < 7. Now, going back to (l47l ). we get Mlz] Q''(/^)] < &c\^^ (N + n)°'^, 
where as stands for the maximum of the left members of the system fl48l) . As < 02, 
we can pick N sufficiently large such that for any n > 0, {N + n)"^""^ > Gc].^, and so 
Ar[z;Q-(/^)] < (Ar + n)"2. 
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It remains to analyze the condition (H). But Ui^^^ = e '^"Vin + {6y)iN^^iN, so that 

^f[mN^^■,c,] <{N+ nr + m <iN+ nr + ci {N + nr-\ 

Now notice that, when m — > oo, we have -, — ptt^- ~ g^^^ remember 

that we have assumed ai — 1 > ^2 — 7, so that if is large enough, the last equivalent 
yields (jv+n+yMA^+n)"! - ^ ^"^^ n > 0. Hence M[yiN^^] Ci] < {N + n + which 
achieves the first step. 

Step 2: contraction property. The contraction argument is now easy to settle. Indeed, if 
^ j^h„ ^(i) — r(^yW)^ then, owing to relation (l46l) . we have 

But U[a{y^^)) - (T(y(2));C0(/^)] < elU[a{y^^^) - a{y^^y)] Q''(/^)], and the previous re- 
lation, together with ([43]) and (SI]), gives 7V[y(^) - y^^^; Q^l^i^)] < c,,^ J^+„7V[y(^) - 

y(2);Q«(/A')],with 

It is finally readily checked that the two conditions 2a2 — (7 — k) < and ai < 0:2 + ! — k 
entail limAr^oo Jn = 0. Therefore, here again, we just have to take sufficiently large for 
the contraction argument to work on the balls S^", n > 0. 

Step 3: patching solutions. The construction of the announced solution y G Qo([0,T]) 
reduces now to a patching argument. Let us make it precise. 

First, define a sequence (y'^,C^")n>o by the recursive condition: (y'',C^') ^ Q'^{Iq) is 
the fixed point of V in and for any n > 1, (y", C^") G Q'^^I^) is the fixed point of T in 

-n-l 

Bn" . This construction is allowed by the first part. Then set, for any t G [0,T], 

Nt Nt 
n=0 n=0 

where Nt stands for the lowest integer such that ^^=0 — ^• 
Ul^_^< s < ... <l§ <t < l^,^^, use the relation 

k'-i 

iSy)ts = e-<'-'^\6y)i., + {5y\i. + e-^*"'" (49) 

i=k 



together with 5x^ = 0, to deduce {6y)ts = a^LCf + f^tsj where r£ = rf/ + rf/ 

fc' 

^tif [''/f '^iN 



' ts ~ -^tl^ 



i=k+l 
k'-l 



i=k 



Owing to the regularity of each C^*, this proves that {y, C^) actually belongs to Qq ([0, T]). 
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Finally, let us go back to the decomposition (14911 to deduce 

fc'-i 

i=k 

Furthermore, invoking the fact that S (^J'{dxz)j = 0, we obtain: 

J,iN {dxa{y)) = J,iN{dxa{y)) + e-^'-^y~^ Jinin {dxa{y)), 

k' —1 k' k' k' —1 

and hence 

k'-2 

{Sy)ts = e-<'-'^^JiN,{dxa{y)) + J,iNjdxa{y)) + Y,e"^'~''-''^Jif^,if{dxa{y)). 

i=k 

Iterating this procedure, we end up with the relation {5y)ts = J^ts{dx a{y)) for all s, t G 
[0, T], which proves that y is a global solution to equation (l22ll . 

□ 

4.5. Application to the Brownian case. We now intend to show that the previous 
results can be applied to a (classical) brownian motion X = {X^^\ . . . ,X^"'^) with values 
in M^'". In other words, we shall consider the processes X^, X^, X^ defined in a natural 
way, according to Remark 14.51 by 

Xl= I e-^'-'^Ux,, Xl = I e-'-'-^'Ux, ® Xl, Xl^ = j e-'^^'^'Ux, ^ {SX^Us, 

J s J s J u 

where X/^ = (pit — v) dX^, and where all the stochastic integrals above are understood 
in the Ito sense. We thus have to prove that those processes satisfy the required regularity 
conditions. 

As far as X^ is concerned, we can use the same proof as in the Young case, and the 
following regularity result is easily shown: 

Lemma 4.17. // is such that /q°° (i^ + < then, for any •y G (1/3,1/2), 

X^ gC7([0,T];R1'") a.s. 

Proof. The same trick as in the proof of Lemma 13 . 1 31 leads to the existence of a continuous 
version of X^{C,) for any fixed ^. Now, just as in the Young case (see the proof of Theorem 
I3.14p . it is readily checked that 6X^ = 0, and thus one is allowed to write, for any p > 0, 

\\Xlm <C\t-s\^ where U,^,,{i) = [[ f^'-^f£', dvdw. (50) 

J Jo<s<t<T \W - V\ 

Our claim is thus easily proved by replacing (1301) with the usual ltd isometry property. 

□ 

Consider now X^, and notice that, up to a Fubini-type theorem, this process can be 
written for all ^ > as 

/■°° . =4 

XLiO= dr^<P{v)X,S^,r^)^XlM, (51) 
^0 

with X^^{^, rj) = e~^^^~'"'> a^u{v)dX^ . The issue then consists in studying the regularity of 

« 4 " 

X . To this end, we shall resort to a GRR-type argument, which requires the introduction 
of a new incremental operator 6 acting on the space C2 of applications on ^2 with values 
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in the space of two- variables functions. This operator should send C2 into the space C3 of 
applications on Ss with values in the space of two- variables functions. 

In order to define 6, observe that for all ^,r7 > 0, 

[aUv) - a™(r/)e-''("-^)] dX, = X,UO««.(^), (52) 

the regularity of which is known. This simple relation yields naturally the following: 

Definition 4.18. If R E C2, let 6R the element of C3 defined by the relation 

{6R)tusi^,v) = iSR)tusi^,v) - atuiORusi^,v) - Rtui^,v)a'usiv) 
for any C,,ri > 0. 

= X 4 

With such a definition, the above relation (l52l ) can be written as: {SX )(^,?7) = 
^tuiO'^usiv)- Furthermore, we have the following equivalent of Proposition 13.121 whose 
proof is postponed to the appendix for sake of readability: 

Proposition 4.19. Let (V, ||.||) a Banach space and fix C,,ri>0. Let R G C2{T]E) such 
that R,X^,v) ^ ^2{V), and set 

J Jo<s<t<T 

where : are strictly increasing functions and 0(0) = 0. Assume now that 

there exists some C{^,ri) > such that, for all ii < £2 & [O;^]? 

'4d(e,r7) 



dtds. 



sup \\{6R)e,uiA^,v)\\<^~' . , 

£i<«<^2 \ \l2 — h 



0(|4-^l|). 



Then, for all s,t E T, 

\\Rtsi^,v)\\<c 



\t^s\ 



^1 Mf7(e,r^) 



d(f){r) 



3. 



It is now possible to give some regularity results for the increment X 
Lemma 4.20. If (p is such that d^\(j){C,)\{l + .^) < 00, then, for any 7 G (1/3,1/2), 

x3 e C3^^([o,r];M' 



a.s. 



= 4 



Proof We will apply of course Proposition 14.191 to X , with ijj{x) = x^^' and (j){x) = 
fjij^g same arguments as in the proof of Lemma [3.131 enable to assert that, a.s., 

= 4 = 

X (^,77) G C2(M '") for all ^,r7 > 0. To find out what C{^,r]) should be, remember that 



s = 4 



{6X )i^uiAi^'n) = ^l^uiO^uhiji)- Hence, according to ((50]), 

||(ll\„,,(e,77)|| < c\l2-u\^{U,,2r,{i)f''W\u-h? 



< c 



1 27 
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Therefore, set 

~ 4 

J Jo<v<w<T \W — V\ ' 

(53) 

and with these notations, 

11^1(^,^)11 <c\t-s\'^ {{U,,2A^,v)Y^''' + {C,,2r,M^,v)Y^'''' 

Going back to (ISB and estimating Xlg{ri) with ||X^s(^7)|| < c\u — s\'^ {U^^2p"iv)y^'^''^" 
for some p" > 0, we get 

II^L(Oil < C\t - uf^ \U - S\'^ R^^2p,2p',2p"{0^ 

where 



R,,2p,2p',2p"{0 = / c?r/|0(r^)|(L/2,,2p'(e,^))'/'^'(f/7,2p"(^))'/'"" 
Jo 

POO _ 

+ / rfr^|0(r^)|(5,,2p,2p'(e,^))'/'^'(t/..2p"(r/))^/^^" 

= -R^,2p',2p"(0 + -^7,2p,2p',2p"(0- (54) 

To prove that G C^'^ a.s, it is now sufficient to show that Af[R^^2p,2p',2p"', ^i] < oo a.s, 
which will be seen as a consequence of E[Af[R^^2p,2p',2p"] ^i]] < oo. 

In order to prove this latter relation, use first succesively Schwarz and Jensen inequal- 
ities to obtain 

E[iU2,M^,v)y^''\U,Mv)y^'''']<Ep^^^^ (55) 

To estimate the first term in the right hand side above, we resort to the fact that 

~ 4 



= 4 _ 

Furthermore, X^^(^, rj) is a random variable in the second chaos of the Brownian motion, 

~ 4 , 

on which all the L^-norms are equivalent. Thus ii^[||X^^(^, r])!^^'] < 1]'^^' \w — vf^ for any 
p' > 1, which yields: 



E[U2y,2p'i^^ V)]'^'"' <CV( II \W- vf^'-'-'^'-' dvdw) 

\J Jo<v<w<T / 



l/2p' 



Hence, if we take p' such that p' — 2'yp' — 1 > 0, that is p' > 1/(1 — 27), then the quantity 
^27,2p'(^5 v) can be bounded as E[U2'y,2p'i^, r/)]^/^^' < C rj. As for the second term of (l55|l . 
we have (remember that n stands for the dimension of X^) 

^[||^l(^)ir]=^ / e-2''("'-"W<n|ti;-^;|, 

J V 

so that the same kind of arguments as for X'^ yield 

E[U^,2p"{v)] <C II k - t^r"-'^*'""' dvdw. 

J Jo<v<w<T 
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By choosing p" > 2/(1 — 27), we get E[U^^2p"{v)] ^ c. Consequently , recalling that is 
defined at equation (l54ll . one gets: 



^[A/'[<2p'.2p"; ^1]] <c 10(01(1 + / \m\v 

Jo Jo 

<c(^^°°rfO0(OI(l + o) • (56) 

As far as -R^,2p,2p',2p" is concerned, we use the definition of C'^,2p,2p', together with the 
previous estimation of E[U^^2p"{v)]i to assert that, if p > 2/(1 — 27), 

i?[C^7,2p,2p'(e,r/)^/^^'f/,,2p"(r/)^/^^"] = r^^i?[f/,,2p(0'/'^^7.2p"(^)'/'^"] 

< r/^E[^,,2,(0]'/'^i?[t/7,2p"(^)]'/'"" <C^^"- 

Hence, 

/■oo /-oo 

EW[Rl2p,2p',2p": C,]] <c di 10(01(1 + / dr^ 

Jo Jo 

<c(^j^ dimm+o) ■ (57) 

Putting together the estimates (l56ll and f l57ll . we end up with E[hf[R.^^2p,2p' ,2p"'i -^^i]] < c>o, 
which ends the proof. 

□ 

It remains to analyze the regularity of X^. To this purpose, we will apply Proposition 
13.121 again, which means that both the moments of X"^ and = X^X^ + have to 
be controlled. We first have to check the following property: 

Lemma 4.21. // drj \4){ri)\ < 00, then, a.s., X^{C) G C2(M"'") for any ^ > 0. 

Proof. This is the same Kolmogorov-type argument as the one used in Lemma 13. 131 The 
details are left to the reader. 

□ 

Let us state now the regularity result for X^: 

Lemma 4.22. If (j) is such that (iO0(OI(l + < then, for any 7 G (1/3,1/2), 
X2 G Cl^{[Q,T]]W'^) a.s. 

Proof. Invoking the fact that 5X^ = X^X^ + X^ and the previous estimations of X^ and 
X\ we deduce \\{5X^),,uiAi)\\ < c\i2 - h?"* mf'^^' , with 



^(^)i/4P5 ^ /?^^2po,2p„2,.(o + (^7,2P3(0)'/'"^ / dr^ \m\{K2vm 



\1/2P4 



We are thus ready to apply Proposition 13.121 to X^ with ip{x) = x^^^ and (j){x) 
^27+i/2p5 to get ||X2 (Oil < c |t - {V'27,4p5(0^^^'" + 1^(0^^^'''}, whcrc we set 



iix^ (nw'^p'- 

' lo<v<w<T \W -V\ 



and so^^[X^;Cl^] < c ^^^^[V2'^Xl■, C,] +Ar[DV4p5.£j|. 
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The fact that M[b^/^P^; Ci] < oo a.s has been shown while studying the regularities of 
and X^, for some well-chosen po,pi,p2,P3,p4. 

To conclude with, let us prove that E[J\f[V2'{X>^,'i ^i]] < C)o: notice that 



Jo 

so that the issue consists in estimating £"[11X^^(^)11^^"^]. To this end, observe that 



E 



4P5 



<c\E 



where X^'i*'^^(0 is defined as e'^^^-^'UX^^ xU^\ But, thanks to the Burkholder-Davis- 
Gundy inequality, we know that 



E[\XS^'\0\'''] < cE 



) \ 2P5' 



< c\w — V 



|2P5-1 



E [{Xl^'^Y'^] ds, 



and 



E[iXl;('^r] = E 



V \J0 



J V 

V 

dr/0(r/)e~''('"*^ dt < 



\s — v\ 



dr](p{ri) 



which gives E[{xI^^^^Yp^] <c\s- and thus £[(X^'i^'^^(0)^^1 < c\w - v\'^'^\ 

In fact, this reasoning remains true for E[{xIv^''^^Y'p^], so that finally -E'[||X^^(^) H^^^j < 
c\w — f 1^^^ If we take p^ > 1/(1/2 — 7), then we get 



E[V2,,,,M)] < 



\w — V 



4p5— 87P5-2 



dwdv < M < oo, 



0<v<w<T 



which leads to the announced claim E[N'[V2^^^^; Ci]] < oo, provided J^^ d^ |0(^)|(1 + ^) < 
oo. 

□ 

We are now able to write Theorem 14. 161 in the Brownian setting: 

Theorem 4.23. Let X = {X^^\ . . . , X^'^'>) a standard Brownian motion on [0, T] with 
values in M}'"^. Introduce coefficients 7 G (1/3,1/2), k G (1/8,7) (^i^d assume that 
/o°^t^^l0(OI(l + O < oo- //^ e C^'\R^''^;W^''^), then, a.s, the system 

'Yo = 

iSY)ts = Jts [dX a (a + d^ Hv) Yiv) 
admits a unique solution in Q''([0, T], R^'"^). 
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5. Application to a fBm with Hurst parameter H g (1/3, 1/2) 

This section is devoted to prove that Hypothesis [3] is fulfilled for a n-dimensional frac- 
tional Brownian motion with Hurst parameter H G (1/3, 1/2). More specifically, we will 
construct a stochastic vector 

(X,X\X^X=^) G L^(fi;C7(M^'") x ^^(Ri'") x ^^^(M"''^) x ^^^(M"'")) 

lying above X (in some rough path sense) such that: 

• X is a fBm with Hurst parameter if, 

.7e(l/3,ii), 

• almost surely, (X^,X^,X^) satisfies Hypothesis [3l that is 

~5X^ = 0, 5X^ = ® + X^ where X^ = 0(0 X^iO- 

Jo 

As mentioned in the introduction, to this end, we shall resort to an approximation of 
the fBm introduced by Unterberger in [23] . Let us recall first briefly the deflnition of this 
approximation in the one-dimensional case. 

All the processes we deal with in the sequel are defined on the same complete probability 
space (n,jF, P). As shown in [23], a simple explicit decomposition of the covariance of 
the fBm allows to introduce an analytic process X+ on the complex half-plane H"*" = 
{x + iy E C : y > 0} such that, if X~ is defined on E" by X~ = Xi , then, for all 

z,w e n+, 



E 



X, X 



w 



E 



X^ Xj7 



0, E 



X^ X- 



w 



2cos(7rif) 



-z{z-w)Y^-\ (58) 



The process X^' has to be interpreted as an analytic approximation of the derivative of 
the fBm, and the simple expression (l58l) for its covariance function is at the core of our 
further calculations. If one desires to construct an approximation of the fBm itself, just 
pick, for t G M and any e > 0, a continuous path 7e^i : [0, 1] such that 7e,t(0) = ie 

and 7£,f(l) = t + ie. Set then X^'^ = ^Xf dz. Likewise, a process Xj~'^ can be 

defined as X^~'^ = J, ^ X~' dz, where ^i;,t ■ [0, 1] H" is such that 7e,t(0) = —ie and 

7e i(l) = t — ie. Of course, X^~'^ = X^'^, and the (real) approximation we shall work with 
is finally defined as 

X^ = 2 Re(X+'^) = X+'^ + X~''. (59) 

The next proposition, borrowed from |23j, gives a first relation between the approximation 
we have just recalled and the usual fBm indexed by M: 

Proposition 5.1. Let X^ be the process constructed above, given by relation ^EM)- For 
all s, t G M, we have 



^rnE[X^X!] = \{\tr + \sr-\t-sr} 



This statement in law will be improved at Theorem 15. 2[ Just notice for the moment 
that, for any fixed e > 0, {X^,t G [0,T]} is a smooth process and (X^)^ = X^-^ + X^_-^. 

In a natural way, the n-dimensional analog of our analytic approximation is a process 
X^ := (X^'(^\ . . . ,X^'*^")), where the components X^'^*^ are constructed from independent 



ROUGH VOLTERRA EQUATIONS 



39 



copies ofX+. We can then introduce the following smooth integral (in the Riemann 
sense) processes associated to X^: 

ft 



ts 



(60) 



J s J u 



dX„= ® {5X 



1,£\ 



The main result of this section, which entails in particular our Theorem 11.11 can be 
summarized as follows: 

Theorem 5.2. Assume that 4> is such that rf^ |0(,^)|(1 + ^^) < oo. Then, for any 
7 e (1/3, if), the sequence of processes (X^, X^'^, X^'^, X^'"^) converges in L^{^1;CJ x C2 x 
^ ^3'^)? ^ tends to 0, to a process (X, X^,X^,X^) . Furthermore, X has the same 
law as a fBm with Hurst parameter H and X^,X^,X^ satisfy Hypothesis 3. 

The proof of this theorem will be carried out in the sections below, and the main step 
in this process will be to prove that (X^, X^'"^, X^'"^, X'^''^) is a Cauchy sequence. This is 
achieved once the following stronger statement is proved: 



E 



M[X' - X"; CJ] + ^^[X^'' - X^-"; CI 



< ce' 



(61) 



with a > 0. To do so, let us fix < r/ < e and follow the same lines as in the Brownian case 
(see Section [475II . which means that the issue mainly consists in estimating the moments 
of any order of the processes at stake. 

For the sake of clarity, the proofs of the lemmas to come are carried over to the appendix. 
Let us also introduce the notation 

X^ = X" - X", X''^ 



X''' - X 



{1,2,3} 



(62) 



5.1. Estimation of the first order integrals. Our approximation results will stem 
from the association of Proposition 13.121 and the following lemma concerning Wiener 
integrals of analytic functions: 

Lemma 5.3. Let < s < t < T and fts an analytic function in a neighbourhood of 

n(s,t) = [z = a + ih eC: a e[s,t\,h e [-{t - s),t- s]} 

such that the restriction fts^i^t] takes value in R. Suppose that fts is bounded on ^{s,t)- 
Then, for any a G (0, 2H), 

21 



E 



fts{u)dXl- / fts{u)dX^^ 



< Cr 



oo>n(3,t) 



\t- s 



,2H-a 



\e — rjl 



where the constant Ca does not depend on s,t,e,r]. 



With this lemma in hand, our approximation result for the first order integrals based 
on X*^ can be written as: 

Proposition 5.4. Let X^ and X^'*^ be the increments defined by (E^. Then there exists 
a constant a > such that 



E 



Af[X' - X"; CJ] + Af[X^^' - X^'"; C 
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for any < rj < e. 

Proof. Recall our notation (i62l) for the differences of increments based on X. For the esti- 
mation of X^, use the classical Garsia-Rumsey-Rodemich inequality to deduce IK^X^)^^!! 

<c\t-s\'' {U^^^pY/^P, where 

which leads to E[M[X^]CT\\ < cE[f/^2p]^/2p_ g^^^ according to Lemma O (take = 
1), we have 

Ernx'^Uf] = E[\\{Xl - X-) - {XI - ] <c\e-r^r\w- v^-^^ , 

for an arbitrary constant ai G (0,2if). Hence, since X+ is a Gaussian process, we 
also obtain (5X^)^^f p] < c |£ - r/l"^^ |u; - so that, if ai G (27, 2if) and 

p> 2/(2/7 -27- E[U[X^]Cl]\ <c\e-T]r''\ 

In order to estimate A/'[X^'^; C^^], notice that we obviously have 5X^'^ = 0, which, 
thanks to Proposition [3l2l gives ||X/j^(OII < c\t - sT {U^,2piOy^^^, where 

JJo<v<w<T \^ ^\ 

Therefore E[Ar[Xi'^; CJ]] < c |0(O|(1+O E[U^^2p{0Y^^'', and invoking again Lemma 

Owith f^,{u) = e'^^"-""), we get E[\\Xl;^{OfP] < c'\e - ri^" \w - Thus, just 

as in the case of X^, we end up with 

POO 

E[Af[X'''^;Cj]]<c\e-vr/' / rf^ 10(01(1 + 

Jo 

which finishes the proof. 

□ 



< c\e — rj 



|"l/2 



5.2. Estimation of the second order integrals. We now proceed to the estimation of 
the increments X^''^ and X^'^, starting with the second one: 

Proposition 5.5. Let X^'^ be the increment defined at ^EMl- Then there exists a constant 
a > such that 



E 

for any < rj < e. 

Proof. As in the Brownian case, write X^'^ as 

xfM) = / dfi 0(/i) x,„ (e, fi) ® xifif,), 

Jo 

with Xj^(^,/i) = e^^^*~''^a„u(/i) dX^. Then the increment X'^'^ defined at equation 
( !62l) satisfies: 

Xttii) = / d^^ 0(^) X,: (e, /x) ® X^f (^) + / d^ 0(/i) X,„ (e, ^) ® X^;^(/i) 

:= (0 (0- (63) 

We will now bound these last two terms separately. 
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4,A 



Let us start by controlling 5^'^'^: notice that {6X = X^''^(^)a(/i), and thus, 

using the same arguments and notations f l53l ) as in the Brownian case, we get 



= 4,A 



= 4,A _ , n-. ~ A,e 



As a consequence, \\X^l ® X^f(/i)|| < c|t - u| |m - 2p,2p',2p"(^' A*)' where 

-^7,2p,2p',2p"(^5 /^) = ^27,2p'(^5 Z^)"^''^^ f^7,2p"(/^)^^^^ + ^7,2p,2p' (^5 /^) ""^^^^ ^7,2p" (A') ^^^^ 

•= -^7,2p',2p"(^5 /^) + -^7,2p,2p',2p"(^5 Z^)) (64) 

which leads to 

^0 Jo 
In order to estimate the first term of the latter sum, start with 

E[ut,2A^,f^y^'''K2Af'y^"''']<EK^^ (65) 

Then use Lemma [5731 to assert that, for any A G (0, 1) and any a G (0,2iJ), 



E 



I re-«(*-")a„s(/i)rfX:- re-«(*-")a„,(/i)dX„^f 



^ 2A I la I, |(2H-q)+2A 

< c/i le — r]\ \t — ' 



Indeed, if 2; G 11(5 j), 



Accordingly, /i)f p'] < c/x^p'Ai |£ - r^]"^^' |w - ^1^'^"°'^'^'^''', and then, if we 

take (Ai, as) G (0, 1) x (0, 2E) such that 2(if + Ai) - aa - 47 > and p' > 2/(2(i7 + 

Ai) - «2 - 47), we get i?[^^^,2p'(e, /i)]'/'^' < c/i^i k - tA'I\ 

To deal with i?[?7^2p"(A*)] (|65l l. consider the following estimation (recall that the 
proofs of all the lemmas in this section are postponed to the appendix): 

Lemma 5.6. Let 0<s<t<T,fi>0. Then 



- l| <2A\z- < 2A \t + i{t - s)- s\^ < cA \t - . 



E 



< c \t — s 



2H 



where the constant c does not depend on s,t,jji^e. 

Therefore, i?[||X^'^(/i)|pP"] < c \w — v\^^^ , which, by taking > — 7), gives 



= A,£ 



-^[^7,2p"(a*)] ^ We can thus assert that -E^i?^ 2p',2p"(^' A^)] ^ \^ ~ V 



7,2p',2p" 



Q2/2 



= A 



As far as i?[ i?^ 2p',2p"(^5 /^)] dMl) is concerned, go back to the definition of C^2p,2p'i 
together with the previous estimations of i?[[/^2p(0] and i?[f/^ 2p"(/^)]) to deduce 

E[X:l,,,A^,^^)] < /^^£^[f/7%(0]'^''^[f>7V(/^)]'^''" < ^^^'\^-vr'"■ 
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Putting together the estimates on and ^R^ we thus have proved that 

< c\\e-7]\'''/^+\e-7]rH. (66) 



Going back to equation fl63l) . let us deal with the term E[J\fpC^'^;Cf^]]. But in this 
latter case, it is readily checked that the previous reasoning and bound fl66l) remain true 

= 4 

by inverting the roles of X and X\ thanks to Lemma [531 and invoking the following 
lemma: 

Lemma 5.7. Let < s < t < T , ^, jj, > 0. Then, for any A G (0, 1), 



E 



\j\-^^'-^^aMdXl\\'' 



< Cfl"^ \t 



2A I . „\2H+2X 



where the constant c does not depend on s,t,^,fi,ri. 

This remark allows us to finally plug the bounds on E[M[^^^^;Cp]], E[^[P<:^^^;Cp]] 
back into equation ( l63ll . and claim that the expected relation 

E[Af[X^'^; Cp]] <c\e- r]f for some a > 

holds true. 

□ 

The upper bound for X^'^ can be written in a similar way as for the previous cases: 

Proposition 5.8. Let X^'^ be the increment defined at Then there exists a constant 

a > such that 

E 

for any < rj < e. 



Proof. Here again, we will proceed as in the Brownian case of Section 14.51 and we shall 
apply Proposition I3.12[ This means that we must control both the regularity of 6X^'^ 
and the moments of X^f . However, ^X^'^ = X^'^ ® X^'^ + X^-" X^'^ + X^'^, so 
that the previous estimations of Propositions 15751 and ISTSl easily lead to || ((5X^'^)t„s(0 II — 
\t - s\^^ D^'^'^iiO, where Z)^'^-'' satisfies E[Af[D^''^''; d]] <c\e- r]\" for some a > 0. We 
have thus obtained that ||Xi^;^(0|| < c |t - {'1^27,4^(0^^^'' + ^"^'"'"(O), with 

J Jo<v<w<T \W - V\ 

and hence E[Af[X''^-C',^]] < c{E[M[{%^^,^)y^P- C^]\ + \e - r^n. 

To study E[^^[{V^^^^pY/^P■,Cl]], we first give a bound on the second moments of the 
increment Xj^j^(^): 

Lemma 5.9. Let 0<s<t<T,^>0. Then, for any A G {0,2H), there exists k > 
such that 

E[\\xf:^{or]<ce^ \t-s\'"-'{i+e), 

where the constant c does not depend on ^, s, t, e, rj. 
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□ 



Since Xfj^(^) is an element of the second chaos associated to the Gaussian process 
X"*", we can easily deduce from the previous lemma that < ce'^^'^{l + 

^2)2^ \w - t,|8^f-2pA_ ^j^^g^ p.^j^ ^ ^ _ ^ ^ l/(4(i7 - 7) - A), 

POD 

Jo 

POD 

Jo 

and the expected result E[Af[X^'^;C2^]] < ce°' holds true. 

We can now conclude this section with the proof of our main theorem. 

Proof of Theorem \5.2. Putting together Propositions l5.4l ISTSl and lSTSl it is readily checked 
that (X^, X^'^, X^'^, X^'^) is a Cauchy sequence in the Banach space L^iVt] Cj XC2 x C^^ x 

The fact that the process {Xf, t G [0, T]} has the same law as a fBm with Hurst param- 
eter if is a direct consequence of Proposition 15. 1[ 

Finally, it is readily checked that the algebraic relations ^X^'"^ = is preserved as e 
tends to 0, by taking L^(n)-limits on both sides of the equality. The same kind of limit 
can be also taken for the relation 6X'^'^ = X^'^ (g) X^'^ + X'^''^, provided one can prove 
that X^'^ and X^'^ are in fact Cauchy sequences in L^(f2). But this is achieved by a 
slight elaboration of Proposition 15.41 thanks to the fact that X^'*^ and X^'*^ are Gaussian 
processes. 

□ 

6. Appendix 

6.1. Proofs of the GRR type propositions. This section gathers the proofs of all 
the general results we need for the regularity of the stochastic processes handled in this 
article. 

Proof of Proposition \3. 1 ^ This is an adaptation of Stroock's proof of the (classical) Garsia- 
Rodemich-Rumsey inequality (see [22j). 

Let s, if: G [0, T] and notice that, for any sequence of decreasing times (sk) G (s,t), 

ms,s,UO = RsAO - Ks,,AO - e-«(^-^'^-+^)i?.,,,.(0, (67) 

so that WRsAOW < \\Rs,^AO\\ + \\Rs,s,^,m + \\i^R)s,s,+AO\\ and by iteration, 

n 

||i?.o.(OII < \\Rs.UO\\ + E + ll(5^)....+is(0ll] • (68) 

k=0 

For all V > s, set I{v) = £ ^/j du, and define the sequence (s^) as follows. 

First, fix Sq G (s, t) arbitrarily. Next, given Sk G (s,t), write Sk = s + Xk {Xk £ (0,t — s)) 
and define < by the relation 2(j){ak) = (f){Xk). Then set Sk+i = s + A^+i, where 
Xk+i G (0, ak) is such that 

I{s + Xk+i) < and i) —— — < . (69) 

«fc \ <p{Xk - Afc+i) / ak 
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Such an element always exists since if we call Ak (resp. Bk) C (0, ak) the set on which 
the first (resp. the second) inequality fails, we have U{^) > J^^ I{s + u)du> ^-^^yu(Afc), 
whereas 



HSk) 



|-Rsfc,s+»(O II 



du> ip 



\Rsk,s+u{0 \\ 

4>{h - u) 



The last two inequalities yield fi{Ak) < ak/2 and fi{Bk) < dk/'^, and thus n^AkUBk) < ak. 
It is then clear that (sk) decreases to s. 

Observe now that 

H^k - A,+i) < HXk) = 2<l){ak) = 4 (^<P{ak) - ^) 



0(A 



< 4 

Plugging this observation into equation (l69ll . we end up with: 

_i f2I{sk 



4(0(afc) - (f){ak+i)). 



< 4 



d(f){r 



< A{(j){ak) - (j){ak+i))i)' 



OikOik-l 



where we have used the fact that ip ^ is an increasing function. Besides, condition (l27l 
entails: 



< 4 



0(Afc) <4V^ 
d(f){r). 



(0(afc) - 0(«fc+i)) 



As -R(0 ^ ^2, we get, by letting n tend to infinity in ( l68ll . 



P.o.(OII <4 / 

^0 

In the same way, we find 

l^*so(OII<4 

Write now 



d(f){r) 



d(j){r) 



to deduce 

11^.(011 < ll^tso(OII + ll^so.(OII + ll(5^)*.o.(OII 



(70) 



< 



#(r) + ||(*fi),„.«)|l, 
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and observe that 

which achieves the proof. 

□ 

We also need to prove a slight extension of the previous proposition to functions indexed 
by two Laplace variables: 

Proof of Proposition \4-19[ It follows the same lines as the proof of Proposition 13.121 Re- 
lation fl67ll has to be replaced with 

which leads to the expected estimation 

whereas (FTOl l becomes 

and thus \\Rts{^,v)\\ < \\Rts,{^,v)\\ + 11^.0.(^,^)11 + \\mtsos{^,v)\\. 

□ 

6.2. Proofs of the complex analysis lemmas. We will prove in this section Lem- 
mas 15.31 15.61 and 15.71 The key ingredients for those proofs are the following elementary 
estimations: 

Lemma 6.1. Let < s < t < T and 7(s^j) the three-part path in 

l{s,t) = [s,s + i{t-s)] U [s + i{t-s),t + i{t-s)] U [t + i{t - s),t]. (71) 
Then, for any e > and any a G (0, 1), 

\dz\ f \dw\\-t{z-w)+e\'"^ <c\t-s\'', (72) 

'7(s,t) "'7(s,t) 

where the constant c does not depend on s,t,e. Moreover, for any rj > and A G {0,a), 
\dz\ [ \dw\U-t{z-w)+e)''-^-i-iiz-w) + ri)''-^\<c\t-s\''-'^\e-ri\^, 



7(s,t) ■^^(s,t) 



(73) 

where the constant c does not depend on s,t,e,r]. 

Proof Denote 7(\^j) = [s, s+i{t-s)], 'y'^^^^-^ = [s+i{t~-s),t+i{t-s)] and 7^^^^^^ = [t+i{t-s),t] 
(see the figure below), so that 



a-2 



/ \dz\ / \dw\\—i{z — w) + e\°' ^= / \dz\ / \dw\ \—i(z — w) + e 

Let us now estimate each term of the latter sum. Notice first that, for any a > 0, the 
integral J^{u + v)"~'^dudv is finite. This allows to obtain: 
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Figure 1. Contours of integration 

Case j = k = 1: 

[ \dz\ [_ \dw\\-i{z-w)+ef-''^ [ [ (u + v + e)"'- 

pt—s pt—s rl rl 

/ / {u + v)''-^dudv<\t-s\" {u + v)'^-'^dudv<c\t-s\'' . 

Jo Jo Jo Jo 



dudv 



< 

Case j — k — 2 



[ \dz\ [_ \dw\ \-i{z -W)+ = ^ /■* - 1^^^ _ + £ _ ^(^ _ ^)|« 

Jf? .^ J J? .^ Jo Jo 



Case j — k — 3 





i>t—s pt—s 


ja-2 _ 






Jo Jo 


< 


\t — s\ . 



^ dudv 



\dz\ I \dw\\—i{z — w) + e\°' ^ 

pt—s pt—s pt—s pt—s 

= / \{t — s — u) + {t — s — v) + 6\"^'^ dudv < / / \u + v\"^'^ dudv 

Jo Jo Jo Jo 

< c\t — s\ . 
Case j — l,k — 2: 







a-2 _ 


ri 




Jo Jo 


< 


\t — s\ 



\dz\ / \dw\\—i{z — w)+e\°' — j j \{u + t — s) — ivl"' dudv 
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Case j = 1, k = 3: 

/ \dz\ / \dw\ \—i{z — w) + e 



0-2 



<(s,t) 



t—s pt—s 



^0 



-i{s — t) + {u — V + t — s)\°' ^ dudv < \t — s 



Case j = 2, k = 3: 



[ \dz\ [_ 



\dz\ / \dw\\—i{z — w)+e 

3 



a~2 



-t—s pt~s 

I \2{t — s) — V — i{s — t + dudv < |t — s|" . 

'0 Jo 

It is clear that the other cases can be dealt with in the same way, which achieves the 
proof of ( l72ll . As for (l73l l. notice that, if z G 7(s,t) and w G 7(s,t), 



- w) + e)""^ - {-i{z -w) + ri) 



a-2\ 



< (2 sup \—i{z — 'w) + x\°' ) (la — 2| sup \—i{z — w) + x\'^ ) \e — ri\ 

x&[ri,e\ x&[ri,e\ 



< c \—i{z — w) + 7] 



(a-A)-2 



SO that the result is a consequence of ([72 



□ 



Remark 6.2. As mentioned in [23], one of the interesting features of the complex analysis 
approach for the stochastic calculus with respect to fBm is that simple deformations of 
contour like (l7T|l allow to transform very singular kernels like (u — into a much more 
tractable term of the form {u + f )°~^. 

We can now begin with the proof of our lemmas, after introducing an additional nota- 
tion: for the sake of conciseness, we shall henceforth denote 



K:{x,y) = i-tix-y) + e^+e2r. 



(74) 



Proof of Lemma 15.31 Write the left-hand-side as 



Ats{e,r]) := E 



n 



t pt 



fts{u)dX^^- / fts{u)dX2\\'' 

J s 

fts{u)fUv) E [((X^'(^))L - {X^'^'X) ((X^'(^)); - (X''-(i));)] dudv. 



Notice then that 

E [(x^'«):(x^'«);] 



E 



{u + le) + X'~'^^\u - le)) [x'+'^^\v + le) + X'-^^^\v - le)) 



= cn{Kl^-\u,v) + Kl 
where ch = T^^^m-, and 

^ 2cos(7riJ) ' 



2H-2i 



E [(X-«):,(X'''«);] = E [(X'''«)UX-«):] = on {K^^-'M + irj^^(^,n)} 
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Using these identities, one can decompose Ats{£,ri) into a sum of terms whose prototype 



is: 



Let us focus then on the estimation of this last term: by a deformation of contour, we get 
(remember that 7(s,t) is defined by f FTTll ) 

Al{e,v) = I dz !dwUz)Uw)[Kl"^-\z,w)-Kl'^-\z,w)l 
hence, owing to ( l73l l. it is easily seen that: 



\A\Ae,r])\ < 



2 

oo,n 



{s,t) 



\dz\ 



\dw\ \K^.'!-Hz.w) - K^^-\z,w) 



^ II xll2 I la u \2H—a 

< c||/IL,n(,,,) k-^l \t-s\ 



Clearly, this argument remains true for the other terms composing Ats{e, if), which achieves 
the proof. 

□ 



Proof of Lemma \5.6l We have 



E 



with E[(X"'«);(X"'W);] = ch {K^^-'^{u,v) + K^l^'\v,u)}, which gives, by an argu- 
ment of symetry. 



E 



ft pt 

c I du dv e-^^'-'^^e-"^'-"^ K^^-\u,v). 



In the latter integral, deform the line [s,t] into 7(^,4) for u and 7(^,4) for v. The result is 
then a consequence of (l72l ). □ 



Proof of Lemma \5. 71 It goes along the same lines as the previous proof, taking into ac- 
count the fact that if z G 7(^,4) U 7(1^, 



\azs{fJ')\ < \e' 



- l|^ < / 1^ - < / |s + i{t - s) - < c/ |t - . 



□ 



6.3. Proof of Lemma 15. 9L As one might expect, the estimation of the variance of the 
convolutional Levy area Xf^'^ (defined by equation ([62l) ) gives rise to more intricate cal- 
culations. We shall distinguish the diagonal and non-diagonal terms, respectively denoted 
by Xfg^'^^'^^ and Xfj^'*-^'^\ and use the fact that 

E[\\x^f{Or]<c{E[\\x^;^'^^^^^ 

The diagonal term X^^^'''"'^'"'^-': For the sake of conciseness, denote := The 
following property of will be useful: 
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Lemma 6.3. Let e,ri > 0, < s,t < T . Then, for any a G (0, 2H), 



E[{X!-X!y]<c^\t-s 



,2H 



X^^f] < c, \e 



where the constants Ci, C2, C3 do not depend on e, rj, s,t,a. 



V\ 



(75) 
(76) 



Proof. The first (resp. second) inequality is a direct consequence of Lemma 15.61 (resp. 
Lemma [53]) • Those two inequalities account for (l76l ). Indeed, we have, on the one hand, 

Em! - X!}' - w - xn'r] < c {ehx! - xir] + eu:^ - x^^y]} <c\t-s\ 

On the other hand, using Holder's inequality. 



iH 



E[{{x! - xn' - {x^ - xn 



2\2i 



Eim 



X?) 



{XI - x^)}' {X! - XI + x,^ - x,n'] 

< E[{{X! - X^,) - {XI - X^s)Y 

and thus 

E[{{Xt - XlY - {X^, - X:ff] < cVE[{X!-X^)^]+E[{Xl-X^r] 

< c\e — ri\ . 

Hence for any A G (0, 1), thanks to (iTTll and (1781 ). we end up with 



(77) 



(78) 



E 



2\2 



A 



= E 
<c\t 



{{XI - Xlf - {X," - X^} 
- s\ |e — ?7| , 



2\2 



E 



l-A 



2n2 



which gives the result if we take A 



Recall from (l60ll that we have Xf/'^^'^^ 



(0 = irdf^Hf^)Xtl ' \^,/i), with 



f 2,e,(l,l) 



X 



ts 



(e,/i) 



due-^^'-'^UXf. / e-^^^'-^'UXf,. 



□ 



(79) 



~ 2, A, (1,1) 

Our main effort will of course concern the estimation of X 



■ts 



Hs 
Hs 



In the absence of 

exponential weights like e~^*^*~"\ namely in the case of the usual Levy area 

dXf^ dX^,, the strategy is obvious (see [23]): one can compute explicitly Af 
|(Xf — X|)^, from which all the useful bounds can be easily deduced. The situation 
is less simple here due to our exponential weights, but we will try to mimic the classical 
situation with a natural trick: integrate the exponential weights by parts. 

More specifically, for a fixed e > 0, since X^ is a smooth process, it is easily derived 
from equation (l79l l that 

Xl'f''\^,^^) = \ {{Xtf - 2e-^^'-^^XIX! + e-«(*-^)(X:)2} 

+ C\e,/i) + ^tf(0 + C(e,/i), 
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where the remainders R , R^' and R are defined by: 

J s J s 

J s 

A little more elementary algebraic manipulations yield then: 

xT^'''\i:^) = \ {{^l - - W - X^^r} + (1 - e-'^^*-)) {XIX! - X^X^} 



+ ke-'^'-'^ - 1) mr - {X^r} + ^(e,/.) := (SO) 



where, of course, R^g{^,fi) = R^^ + Rts (0 + Rts (^^A^)- Recall that we have to 

estimate and according to (l80l l. we will treat the diff'erent terms Kj 

separately. 

Study of Ki : The expected value E[Kf] can be bounded easily thanks to ( 1761 ). applied 
for some fixed a G (0, 2H). 

Study of K3 : By Holder's inequality, we have 



E 



< e\t-sfE[{XI-X^)'{XI + X^y] 



< e \t-s\' Vmixi - x^rWEiixi - x^nxi + r^y) 

< ce\t-s\^\e-r,f\ 

where we have used (l75l ) to get the last inequality. 

Study ofK2 : Likewise, using the decomposition X.^Xf-X.^X^'' = (Xf-Xt'')X,"+Xt''(X,"- 
X^), we get, for any /ii,/i2 > 0, 

E[{1 - e-^^(*-"))(l - e'^"^^'-'^) {X^X! - X^X^}] 

< f,,f,, \t-sf {E[{X! - X^)Xt] + E[Xnx: - Xm (81) 

< \t - sf{y/E[{X! - x^y] + ^E[{xi - x3r]] 

= A = A 

Study of K4 : For our purposes, it remains in fact to estimate E[R^^{^ , ^i) R^g{^ , ^2)] for 
any /ii,/i2 > 0. This expression can be decomposed into 

ijG{l,3} fcg{l,3} ie{l,2} 
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whose terms will be treated again separately. For the last term, we have, just as above. 



< 

- 4 



< ce\t-sf\e-r]r/\ 
By the same arguments, we easily deduce 

E[Rtf{ORts\^, /i.)] <ca^-f^Mt-sf\e- vr^' for ^ e 1, 2. 

To deal with E[R^g (^,/ii)_R^^ (^,/i2)], use the same trick as in (ISTll to deduce 

E[lCf{^,fii)Rt'\^,fi2)] 

<ie-/^i| 1^-/^21 [ du [ dvE[{xtx^^-x^x:i}{x^,x^,-x^x:!}] 

J s J s 

<c\i- /Uil 1^ - /i2| \t - \e - rf''^ . 

« 1,A 

In order to handle the terms involving R^^ , we resort again to the integration by parts 
method, which yields: 



where 



Rts = c|/ili,(^,/i) + /i(^-^)Si,(^,/i)|, 

Xs{i.i^) = -^l t due-^^'-^^Xl^ \\ue-^^'~-\Xlf 

J s J s 

b]m, /.) = du e-«(*-")x: ^" dv e-'^^^-^^X.^ 



Then it is readily checked, by some elementary computations, that: 

2 

E[Rt:\i,^l^)IC:\i.^^^)] < c\t-sf\e-vr/'{i+e}ii{i+f^i} 

k=l 

E[Rt:\^,fii)Rt'\o] < c\t-s\'\E-vr^'{fi,^+fi.e+f^u} 

2 

1=1 

It is also easily seen, by means of the same considerations, that: 

2 

E[XT'^'''\^, /^2)] <c\t- .|^^-° k - {1 + n {1 + /^'} ' 

1=1 

for some k > 0. 

Finally, gathering all the estimates we have provided so far, and assuming the condition 
dfi |0(/i)| {1 + /i^} < oo on 0, we end up with the announced estimation, that is 

^[(^2A(l,l)(^))2] < ^ 1^ _ ^^4H-a^. {l+e}, (82) 

for a certain strictly positive k. 
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The ofF-diagonal term X^'*'"'^'^'': Recall that Xj^j*-^'^-* is defined by: 



Js 



Hence, the moment £'[(Xi'^'*'^'^-*)^] can be written as 

oo POO Pt pt PXl PX2 

dyi<P{yi) / ^2/2 0(1/2) / / / / ^3;idx2rfx3c/x4M^;^2X3X4' (83) 

Jo JsJsJsJs 

with 

(^''•^^^)k(^'''^'^)L3}{(^^'^'^);.(^^'^'^);.-(^^'''^'^)L.(^'''^'^);j]. 
But according to our convention (l74l ) for the fractional Brownian kernel, we can write 

E[(X^'«);^(X^'(^));3(X^'«);,(X^'(^));j (84) 

= ^ [(^^•^^^)k(^^'^^^);j ^ [(^^•^'^);3(^^'^'^)^J 

= 4 {irJf-2(xi,X2) + irjr'(x2,a;i)} {irJf-2(x3,X4) + irJf-'(x4,X3)} , 

where we have used (l58ll in order to compute expressions like -E'[(X^''^^))^.^(X^''^^))^.J. In 
the same way, one can check that 

^[(^e.(l)),^(^.,(2)),^(^.,(l)),^(^.,(2)),J 

= 4 [Kl'^-^x.^x^) + Kl'^~\x2,x^)] {irJJ-2(x3,X4) + Kl^-\x,,x,)] , (85) 

so that the difference between Expressions (l84l) and (l85l) . denoted by A^'^''^^^^^^, can be 
decomposed as A^-'' = c]j{A'''> + 5^''' + C'''> + D^'''}, with 

^^f...3.. = {^Jf-'(^i'^2) - Kj?-^(xi,x,)} Kl^-\x,,x,) 

Bt;^x2x,x, = {Kl'',-\xuX2)-Kl^-\x,,X2)]Kl'',-\x,,x,) 

+ -^e,^ ^(2^1,2:2) ^(3^4, 2:3) — ^(x4,X3)} = B^:^^^^^^^ + Bl;^^^^^^^, 

= {Kl^e-\^2.x,)-Kl^-\x2,x,)]Kl^,-\x,,x,) 

+ irjj-2(x2,xi) {K2f-2(^^^^^) _ Kl^-2{x:,,x^)] , 

^^f™ = {i^'.'.r'(^2,xi) - K^;^-\x2,x,)} K^^-\x„Xs) 

+ K^^-\x2,x,) {K^j-\x,,xs) - K^^-\x,,X3)} . 

Thanks to Lemma [6?n the treatment of A'^'^ in the expression (l83l) becomes easy. Indeed, 
we are allowed to deform the expression 

A'''''^= [ dxi [ dx2 [ dxs [ c/x4e-«(*-"^)e-^^("^-"^)e-«(*-"^)e-^^("^-"*)i^f43,^ 
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into 



I dz2 I dw2e-'''^''~'^^e-''^^'"'-'"'^K^^~\z2,W2), (86) 

where 'j(s,t){z) stands for the path ^(s^t) stopped at z, defined similarly to Proof of 
Theorem 3.4]. Hence, invoking (l72l l and ( TTSl l. we get 

< / \dz\ [ \dw\\K^J~\z,w)-K^_^-\z,w)\ 

I \dz\ I \dw\\Kl^,~\z,w)\<c\t-st"-''\e-r^\'' . 

The same arguments hold for A"^'^'^, as well as for D^'"^. 

The estimation for B'^'^ and C^'^ is less obvious, since we must cope with integrals of 
the form 

B^'^'^ = fdxi rrfx2e-«(*-^i)e-«(*-^^) {K'^^-\x,, X2) - K^^-\xi,X2)} 

J s J s 



X4,X3j, 

for which the complex deformation (l86l) is not allowed {K^^~'^{w2, 22) would be ill-defined 
for small e, since —i{w2 — Z2) might be negative, see |23] for a further explanation). In 
fact, the result is a consequence of the technical lemma below. Indeed, with our notations, 
( l88ll is equivalent to 

<c\t- sf'- e-{l + e){l + + /i^), 

which, as Jq°° ci/i |0(/i)|(l + /i^) < 00, gives an accurate bound for our purposes. To 
conclude with, it only remains to observe that the reasoning which leads to (i88l ) can be 
easily adapted to 5^''''^. The term (7^''^ is then handled with an argument of symmetry. 
Putting all our estimates together, we have thus proved that 

E[{Xl^^^'^'\Of] <c\t- sr-e- {l + e). (87) 

for a certain a > 0. 

Owing to inequalities ( l82ll and (l87ll . the proof of Lemma [5^91 is now easily finished. We 
are thus only left with the proof of the following lemma: 

Lemma 6.4. Let s <t, ^,/ii,/i2 > 0, and set 
Qif (e,/ii,/i2) = fdx, rf/x2e~«(*-^'^)e-«(*-^'^) {K',;^-\x,, X2) - K^^-\xuX2)} 

J s J s 



Xi l>X2 



X4,X3j. 



Then, for any a < 4H — 1, we have 

2 

|Q?f(e,/^i,/i2)| <c|t-s|^^-"£"(l + an(l+/^')' (88) 
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where the constant c does not depend on s,t, e^rj, ^, fii, fj,2- 

Proof. First, notice that the estimation is obvious if t — s < 2e, since then \K^^~'^{x4,X3) 



< \t — ^ and 



\T^2H-2/- \ T^2H-2/ M ^ 2H—2—a\ I" ^ Ij. \2H-2-a a 

\K^^ ''{xi,X2) - K^^ {xi,X2)\ <cs "\e-ri\ < c\t - s\ e . 

From now on, we thus assume that 2e < t—s. The strategy in order to control our multiple 
integral Ql's{i, jJ^i, /U2) consists then in two main steps: (i) Handle the exponential weights 
by means of successive integrations by parts (recall that is a smooth process for a given 
£ > 0). (ii) Control the singularities of the fBm kernel by a convenient deformation of 
contour. 

A first application of the integration by parts trick gives the following identity: 

dx3 / dx4e-^''^^'-'''^e-''^^"'^K^^-^{x4,Xs) 

J S 

= ciji^jf ^^^^ ^^^) _ e-'^i(^i-^)ir2J(a;2, s) - e-'*2(x2-s)^2^ ^.^^ 

where R can be further decomposed into Rxix2{i^1j 1^2) — -^iiX2(A*i' A^s) + -Rxix2(a*i) + 

rxi px2 

dXs J dX4 e-^i(^i-^3)g-M2(x2-X4)^2^-l(^^^ 

dxse-^'^'-''''^K^,^{x2,X3) 
iHi J dxs e-'^i(^i-^3)g-M2(x2-s)^2H 
We have thus proved that 

Qt'J'it /^l, 1^2) = ItsiO + ^hsit /^l, 1^2) + llltsit Hi, H2) + l^tsit 1^1: /^2) + C2-RxiX2(/^l, /^2), 

with 

4(0 = J^xi jyx2ml'l[,^{0K'^;{x2,xi) 
I%,(e,/xi,/X2) = [ dxi f dx2m'^l^^{i)Bl\^^{ix^, 112) 

J S J S 



^xiX2(/^i'/^2), with: 




^ilX2(/^l'/^2) 


= C//2 


-^XlX2(/^l) 




^IlX2(/^l'/^2) 


= 



int,(e,/ii,//2) = ^ ^^1^ dx2rrf^l^^{i)Cll^^{ni, H2) 

IVts(^,/^i,/^2) = j^^^^ dx2ml'l'^^^{^)Dll^^{ni, H2), 

where we have set /^^'^^^^(O = e-«(*-^i)e-«(*-^2)|/^2^-2(^3,^^ 3,^^ _ ;;^2H-2(^^^^ 

now estimate these 4 terms separately. 

To begin with, let us consider the case of ItsiO- elementary change of variables 
{u — Xi — X2, V — X2) yields: 

4(0 = r ' dufliu) {{-tu + 2£)2^-2 - {-zu + e + rjf''-'} {zu + 2^)2^, 

J-{t-s) 
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where 

fi{u) = l|„<o} / ^it;e-«(*-"-^)e-«(*-'') + l|„>o} / dv e-^^'^^-'^^ e-^^'~'\ 

J s~u J s 

It is easily seen that for u E (0, t — s), 
so that 

4(0 = 2 Re (^j^ ' du fl{u) {{-tu + 28^"-^ - {-lu + e + r^)^^-^} (lu + 2^)=^^ 
Write the latter integral as 

J[-2ie,i-s-2ie] 

and deform the line \—2ie^ t — s — 2ie] into a four-part contour ^^{t-s) = 7^ (^t-s) ^ 7e (t-s) ^ 
Is (t-s) '-'^e (t-s)' '^here 7^ runs along the half-circle centered at the origin from —2ie to 
2ie in {z : Re 2; > 0}, 'j'^ (t-s) ^^^^ [2^£^, '^(i — s)], 1^ (t-s) li'^^ [i{t — s),t — s+i{t — s)] 

and 7s the line [t — s + i{t — s) ,t — s — 2ie\. 



i{t - s) 



2ie 



r 



t - s 



-2ie 



Figure 2. Deformation of [-2ie, t - s - 2ie] 

Using the decomposition 

e-«« - e~i('-'-^)e-^^'-'^ = e-^(*-"-")(l - e~^^'"^) - (e-^^*"^"") - 1) + (e'^" - 1), 

it is readily checked that sup^g^^ Iftsi^ + 2i£^)| < c |t — s|. The estimation of (iQOll on 
each of the 7* is then a matter of elementary calculations, that we proceed to detail 
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now: for 7^ (t-s)^ if 2/ ^ Iv^ ^] 9 G [— vr/2, 7r/2], that | — 2ie e*^ + 3e + y\ > e. 

Thanks to the fact that e < (t — s)/2, this leads to 



[ dz fi{z + 2ie) {{-iz + 4^)2^-2 _ ^ 3^ _^ ^)2//-2| .2// 



< C|t-S|£^^"^ <c\t- S 



AH -a 



For the upper bound on the path 7^(j_g), notice that if x G ^e,t — s], ?/ G [r/, e], then of 
course |x + + ?/| > x. Thus, for a small enough positive parameter a, we obtain: 



I dz fi{z + 2ie) {{-iz + 46)2^-2 - {-iz + 3£ + r^)^^-^} 



2H-21 ^2H 



t-s 



2H /.2H-2-a 



< c\t — s\ dxx {x 

J2£ 

< dt-sl'^-^e" C dxx'""-^-^ 



< c\t — s\ e 



since, by hypothesis, 2 + a — AH < 1. For 7^ start with |t — s — m + Se: + ?/| > |t — s| 
if M G [0, t — s] , ?/ G [r], e] , to deduce 

j dzfi{z + 2ie) {{-iz + 4^)2^-2 _ ^ 2^ + 77)'^^'} 

ft — S 

< c\t — s\ du \t — sf^^'^^°' \e — ril'^ \i(t — s) + uf^ < c\t — s\'^^~°' e"". 



Finally, as far as 7^(t_g) is concerned, observe that for any v E [0,t — s + 2e], y E [i],e] 
|— i(t — s) + (t — s) — V + 3e + y\ > \t — s\, and thus 



%-(t~s) 



dz fl{z + 2ie) {{-iz + Ae) 

t-s+2e 



2H-2 



-IZ 



z 



< c\t-s\ 

Jo 

< c\t — s\ e . 



dv \t — sf^ ^ ° |e — ?7|" \t — s + i{t — s) — iv 



2H 



Therefore, these four elementary bounds, computed on the paths 7^(-j_g) up to '^'^ 
allow to claim that \its{i)\ c\t — s\^^~'^ e°-. 

Consider now the term lYts{i, 1^2)'- it is readily checked that 



t rt 



- X2) + 2ef^-^ - {-i{xi - X2) + 5 + rif^'^] 
and perform the same deformation as in Lemma [6. II to deduce 
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In order to deal with llts{^, /J^i, IJ>2), notice that 



and write, by means of another integration by parts, 

y"*cia;ie-«(*--)e-'^^(--) {Kf,-\x,,x,) - Kl^,-\xux,)] 

= 4 



= c 



with 

=:4 



J s 

Now, use the same strategy as for lts{^-, Hi, which consists here in writing 
^* dx, e-'^^(*-)e-^(*--)X,2J(a;2, .) x,) - K^^^-\t, x,)} 



— e 



-fii{t-s) 



' [2ie,t~s+2ie] 

{{-i{t -s-z)+ Aef"-^ - {-i{t -s-z)+3e + 
and deforming the line [2is, t — s + 2i£] into 

%,(t-s) = ll,(t-s) U [-Sie, -i{t-s)]U[-i{t-s),t-s-i{t-s)]U[t-s-i{t-s),t-s + 2ie] , 
with 7^ the half-circle centered at the origin from 2ie to —2ie in {2; : Re^; > 0}. The 
same kind of elementary estimations as for Ifs(^, /Xi, then lead to 

Jjx,e-^^^'-^^e-^^'-^^^K^;{x,,s) {K^^^'it^x,) - K^,^,-\t, x,)} 



< c\t — s\ e . 



Likewise, 



/ dx2 e-«(*-)e-«(*--)ii',^,f (0:2, .) X,) - K^"^-\s, x,)} 



and accordingly 

Those arguments can be easily adapted to IIIts(^, A*i, A*2) to get 
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W 



.'here 



We finally have to cope with the remainders R, R ,R . Owing to the higher regularity 
of those terms (as regards the kernels), it is rather clear that simple integration by parts 
should be sufficient to reach the expected bound. 

^ 2 

Consider for instance the case of it! ni) defined by: 

= c/xi ^* dx dy e-^(*-)e-«(*-^) {i^s,f "'(^, v) ' v)} v)^ 

with 4>^^{x,y) = due~'^^^^~"^K'^^(y,u). Another integration by parts yields: 

dy e-«(*-^) (x, y) - Kl^-\x, y) } ct>'^^ (x, y) 

{Kl^,-\x,t)-Kl^^-\x,t)]rKx,t) 

{Kl^-\x,s)-Kl^-\x,s)]r{x,s)\ 



-c j dy{K'^^r\^,y)-K',!i-\^,y)}e-^^'-'^\^ir'i^,y) 

~ 2 

and thus, plugging this expression into the definition of it! and integrating by parts again, 
we obtain 

RiCf^i) (91) 
^c^,^[[{ {2ef" - (£ + nf"] {t, t) - e-^(*-) {K'J (s, t) - K^^ (s, t) } {s, t)] 

- e-^^*-^) [ {i^J?(t, s) - Kl^it, s)] r\t, s) - e-«*-) {(2^)^^ - {e + r,r] 0-(., .)] } 

+ R'\i,^xl) + R'\i,^x^), 



with 



t rt 

dii n 

xy 



R'\^,fx,) = c/XiM dx I dyql'^^{Or'{x,y) 



dy 



R {^,l^i) = I dx I dyql.'^{^)^^{x,y), 



s •/ s 
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where we have set q^yiO = ^ 



{K^^-^{x,y) - K^^-^{x,y)}. Let us now 



integrate those last two expressions by parts with respect to x: we obtain 



= 2,1 



=ci/ii^ / dye 



-at-y) 



{K',;{t,y)-K^^^{t,y)}r'{t,y) 

- e-«(*-) {irjf is, y) - Kl^ is, y) ] is, y) 
+ C2 /ii e f dx I' dy e-«(*-)e-«(*-^) {Kjf (x, - K^^ (x, } (x, y) 
+ C3/xie Trfa; f dye-^^'-^^e'^'^'-y^ {Kl^,ix,y)-Kl^,ix,y)]^ix,y), 



s 2,2 



and the expression for R is 



=ci/ii / dye 



AMI 



{Kl^it,y)-Ki:iit,y)]^-^it,y) 



e-^^^-^){Kl^is,y)-Kl^,is,y)]--^is,y) 



+ C2/iie/ dx I dye-^^'-^^e-^^'-y^{Kl^,ix,y)-Kl^ix,y)]^ix,y) 



+ C3/ii/ rfx/ dye-«(*-^)e-«(*-^){irJf(x,y)-irJJ(x,y)} 



dydx 



ix,y)- 



The following (easy) estimations come then into play: whenever x,y E [s,t], we have 



AMI 



2^+1 



AMI 



(9x 



< c |t — s|^^ {1 + /ii} 



AMI 



5y 



< C |t — sf^ {1 + yUl} 



and if a G M, \K^^ia,0) - K^^ia,0)\ < ce"^""^ \e - r]\ < ce"^^. Besides, 



dydx 



[X. 



y) = c |Atr^(y, x) - /ii du e-^^(--")irJf-^(y, «)| . 



Going back to (I9T1) . the previous estimations finally give rise to 



= 2 



which leads to the expected bound since 2H + 1 > AH — a and 2H > a. 



= 4 =5 



The same arguments enable to handle R ,R ,R ,R , which achieves the proof. 



□ 
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